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Chapter 1 

DISCUSSION OF SIGNIFICANCE AND SCOPE 
OF THE UNDERTAKEN STUDY 


1 .1 Introduction 

Research and development in advanced composites technology and 
design procedures has advanced to a point where these materials can 
now meet the challenge of highly demanding aerospace, automobile and 
other structural applications. Performance payoffs in the form of 
weight savings, component integration and a five-fold maintenance re- 
duction as compared to metal design have been demonstrated (Forsch 
[1]) and contribute to the overall picture of an attractive marketable 
product. 

However, resin matrix composites have moduli and strength 
properties which are viscoelastic or time dependent. As a result there 
is a need for methods by which short term test results can be used to 
predict long term results. 

The current effort is part of a continuing cooperative NASA-Ames - 
VPI&SU research program directed towards the development of an 
accelerated strength and stiffness characterization procedure which 
generates the design data and the confidence levels, needed to insure 
structural integrity for any projected lifetime. This process is 
greatly complicated by a number of environment exposure variables 
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such as temperature, moisture*, UV and space radiation, high vacuum 
and many others. An examp l'e of a typical thermal profile for an air 
carrier given by Ripley [2] is given in Fig. 1.1. Any attempt to 
handle such variables by overdesign will penalize the economical 
features mentioned above as a knockdown factor will seriously limit 
the full weight reduction potential of advanced composites. 

1 .2 Previous Efforts 

An accelerated characterization, which is schematically given in 
Fig. 1.2, has been developed and documented by Brinson et al . [4-13]. 

A fair amount of experimental support has been accumulated, based on 
elements of linear viscoelasticity by Yeow [5]. Nonlinear visco- 
elasticity came into this picture when Griffith [7] potentially identi- 
fied this behavior as the reason for a serious underprediction of 
creep rupture data for [90/±60/90]2 $ laminates. The prediction was 
based on a linear incremental lamination theory developed by Yeow [6]. 
As a result, a study of nonlinear viscoelastic characterization pro- 
cedures was performed by Griffith [7] which led towards a valid 
graphical Time-Temperature-Stress Superposition (TTSSP) principle for 
the generation of modulus master curves. 

Dillard [9] developed a non-linear viscoelastic incremental 
lamination theory and found reasonable correlation between his predic- 
tions and experiments. Because of the limited experimental data 


*It has been shown that moisture pickup is generally reversible. 
When combined with thermal spikes, however, this phenomenon can create 
irreversible damage as shown by Adamson [3]. 
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Fig. 1.2 Accelerated Characterization Method for 
Laminated Composite Materials. 
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available from the work of Griffith, Dillard used a semiempirical non- 
linear viscoelastic model due to Findley to represent lamina 
compliances. 

1 .3 Objectives and Outline of the Present Study 

The objective of the current study was to cast the graphical 
TTSSP method into an analytical form. It was believed that the 
analytical framework developed by Schapery in a series of publications 
could be adapted for this purpose. A literature review on nonlinear 
viscoelastic constitutive equations is given in Chapter 2, while the 
thermodynamic basis and the derivation of the Schapery model are the 
subjects of Chapter 3. 

An effort to relate viscoelastic deformation and failure on an 
energy basis is developed in Chapter 4. This allows an immediate trans- 
fer of any modulus master curve into a strength master curve as 
depicted in Fig. 1 .2. 

Chapter 5 concentrates solely on 934 neat resin properties and 
on the discussion of obtained experimental results thereto. 

Master curves for transverse and shear properties of T300/934 
as obtained through the Schapery analysis are discussed in Chapter 6. 

Results on creep rupture for [90]g s and [60 ]q s lamina, based on 
the energy criterion (Chapter 4), together with a comparison of an 
incremental lamination theory prediction with experimental results on 
[90/±45/90]2 s are discussed in Chapter 7. 

A general discussion and conclusions are given in Chapter 8. 
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Chapter 2 

INELASTIC TIME-DEPENDENT MATERIAL BEHAVIOR 
2.1 Introduction 

The purpose of this chapter is to produce a brief literature 
overview on research in nonlinear viscoelasticity. The majority of 
publications explain how to deal with magnitude nonlinearities , i.e. 
they account for the stress deoendence observed on creep compliance 
versus time as illustrated in Fig. 2.1. 

A much smaller number of articles also account for i ntermode 
non! inearity (nonlinear interaction between different components of the 
stress tensor) and almost invariably use a single mechanical in- 
variant, namely the octahedral stress. In this respect it should be 
mentioned that data generated by Cole and Pipes [14] shown in Fig. 2.2 
indicate the extent of stress interaction in off-axis unidirectional 
boron-epoxy coupons. These experimental results show that inelastic 
behavior is not a monotonic function of fiber orientation. Similar 
behavior has been observed for graphi te-polyimide off-axis coupons by 
Pindera and Herakovich [15]. They indicate that the inclusion of 
residual thermal stresses and axial stress interaction at the micro- 
mechanics level are necessary in order to predict the correct trends 
in the shear response of off-axis specimens. 

An almost never mentioned type of nonlinearity is i nteraction 
nonl inearity , i.e. nonlinear interaction between events (i.e. load 
application or release) occurring at different times. 
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2.2 Representations of Constitutive Equations in Nonlinear 
Viscoelasticity 


2.2.1 Multiple Integral Representations 

The multiple integral representation is a very appealing 
theoretical concept, since it is not limited to a particular material 
or class of materials. The Volterra-Frechet expansion for a one- 
dimensional case is written as: 
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't 

0 ^ (t - x i )da(x-j ) + 

—00 




T 1 




x 2 ) da ( x i ) da ( x 2 ) 


+ . . . 


( 2 . 1 ) 


The first integral of Eq. (2.1) describes the linear viscoelastic 
behavior, defined by the Boltzmann theory and the second and higher 
order integrals are representations of both magnitude nonlinearity and 
interaction nonlinearity , the latter implying an interaction effect 
between, e.g., the stress increments at times x^ and X2* 

If a force is suddenly applied to a specimen, producing a constant 
stress a, it follows from 2.1 that 

6 = D 1 (t)a + D 2 (t,t)a 2 + D 3 (t,t,t)a 3 + ... (2.2) 

The functions D^(t,t,. . . ,t) can thus be determined from creep tests. 

The number of stress levels must be equal to the number of kernels J^. 

Gradowczyk [16] showed, however, that the system of linear equa- 
tions which serves to determine the unknown functions is ill- 
conditioned and the procedure of finding the kernels is unstable. 
Being unaware of Gradowczyk 1 s work, this author attempted to use the 
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multiple integral representation for p.v.c. and had to give up the 
attempts due to the extreme instability of the results obtained [19]. 

The extension of the Volterra-Frechet expansion to the 3- 
dimensional case has been given by Green, Rivlin and co-workers [17]. 
Even if we limit their expression to triple integrals (small non- 
linearities) twelve kernels still must be determined by experiment. 

An alternative of the Green-Rivlin approach has been formulated 
by Pipkin and Ropers [18] and is often referred to as nonlinear super- 
position theory (NLST) which can be written as 


e 



D-| (t - 


t , a(x ))a(x )dx 



t-| ,a(x 1 ) ,t 


- X 2 ,a(x 2 ) )a(x-| )a(x 2 )dx^ dx 2 + ... (2.3) 

It is clear that eq. 2.3 is nonlinear even on the first approximation. 
Thus, the additivity of incremental stress effects in the Boltzmann 
superposition sense is preserved. 

Kinder and Sternstein [20] showed that difference in strains pre- 
dicted by the two theories is in the transient portion if the material 
exhibits fading memory. Thus, in a single creep test, the equation 2.3 
reduces to 


e(t) 


; D,(t - T, 3 (-))5(-)d- 


(2.4) 


which is the most general representation of nonlinear creep provided 
the transient response from the step loadings has decayed. However, 
whether multiple loading and/or unloading programs require additional 
integrals to describe time interactions of the loading history remains 
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an open question. 


2.2.2 Single Integral Representations 

i) Modifications of the Boltzman Theory of Linear Heredity 

Two schools of thought have been identified in the literature. 

The first one states that it is, in principle, possible to map a set 
of nonlinear isochronous stress-strain curves into a linear set, pro- 
vided we nonlinearize the stress and strain measures as 

»[«(»)] ■ ( D(t - t) jr .[-(-) ]d- (2.5) 

* “ Cfj 

Here the equation is linear but only in the nonlinear stress and strain 
measures j and t. This form has been used by Koltunov [21]. A simpli- 
fied form of Eq. 2.5 have been proposed by Leaderman [22] 

6(t) = ; D(t - t) ^ .[-(T)]d- (2.6) 

' - oo 

and by Rabotnov [23] 

»[e(t)] = ' o(t - t) Si d- (2.7) 

J -CC 

The second school allows for a stress or pressure-deformabl e time 
measure. Indeed, the pressure dependence of the mechanical response has 
been predicted theoreti cal ly on the basis of the free volume approach 
by Ferry and Stratton [24] who derived the pressure analog of the well- 
known WLF equation in the form: 



where a„ is the ratio of the relaxation times at pressure p to the 
P 

relaxation times at pressure p Q , B is a constant assumed to be of the 
order of unity, f is the fractional free volume at the reference 
pressure, and is the isothermal compressibility of the free volume. 
Beuche [25] and O’Reilly [26] proposed the equation 

log a p = C(p - p Q ) (2.9) 

where C is a constant and other variables are as defined above. 

Knauss and Emri [27] postulate a linear dependence of fractional 
free volume, on the instantaneous values of temperature, T; water 
concentration, C; and mechanical dilatation e, i.e. 


f = f + aAT + yC + 66 (2.10) 

where a is the coefficient of volume expansion, y is the moisture co- 
efficient of volume expansion and 6 is the coefficient of mechanical 
volume expansion. By substitution of (2.10) in the Doolittle equation 
[28] 


log a = 


B 


2.303 



( 2 . 11 ) 


they obtain the following very general temperature, moisture and pressure 
dependent shift-factor 


log a = 


B 

' 2.303 f 

o 


aAT + yC + 69 

f + aAT + yC + 69 


( 2 . 12 ) 


Neither equation (2.8), (2.9) or (2.12) has so far been subject to 
extensive mechanical testing on solid polymers. 
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ii) Single Integral Representations Based on Thermodynamics 

In order not to obscure the main issue of a literature review, 
we shall not attempt a full treatment of the thermodynamic foundations, 
but rather concentrate on the hierarchy of different continuum thermo- 
dynamic theories. This hierarchy, shown in Fig. 2.3, has been reviewed 
in a wel 1 -documented survey paper by K. Hutter [29]. Figure 2.3 shows 
that the main concepts which can be utilized to develop constitutive 
equations for inelastic materials are the rational thermodynamics 
approach and the state variable approach (also called irreversible 
thermodynamics). The latter will be discussed in detail in Chapter 3. 
Coleman and Noll [30] applied functional analysis to continuum mechanics 
in order to express current stress, free energy and entropy 
as functionals of the thermo kinematic history. Their analysis 
continues to influence rational thermodynamics research up to this date. 

The state variable approach includes certain internal variables 
in order to represent the internal state of the material. Constitutive 
equations which describe the evolution of the internal state are in- 
cluded as part of the theory. Onsager [31] published the concept of 
internal variables in thermodynamics. This formalism was used by 
Biot [32] in the derivation of constitutive equations of linear visco- 
elasticity. Ziegler [33] extended this work to include nonlinear 
viscoelastic materials and his approach has been used by Schapery [34] 
and Valanis [35] to develop specific theories. The extension to in- 
clude viscoplastic materials was made by Perzyna [36]. 
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2.2.3 Mechanical (Rheological) Models 

Elastic behavior is compared to that of a massless Hookean spring 
and the viscous behavior to that of a Newtonian dashpot. The success- 
ful description of the stress strain behavior depends on the proper 
choice and arrangement of the elements or hypothetical rheological 
units. The final hardware-system that describes the behavior of a 
material system in the linear range can, at least in principle, be 
extended into the nonlinear range by the introduction of nonlinear 
springs and dashpots. An example is the nonl inearized generalized 
Kelvin element whose mathematical equation for creep can be written 
as, 

„ f 2 (o) N 1 

6 = ? f-i(a) + — — — • t + fo(a) I f- (1 - exp (-t/x.)) (2.13) 

L o 1 n i=l L i 1 

where f. (a) is a nonlinearizing function of stress. Bach [37] used 
this form in his studies on wood. He separated the total creep strain 
into instantaneous, delayed-elastic and flow components. The 
instantaneous or glassy response was always independent of stress 
(thus f-| (a) = 1), but the delayed elastic and flow compliances were 
independent of stress only within certain limits of stress, depending 
on moisture content and temperature. 

We will show (Chapter 3) that Schapery's thermodynamic model, at 
least under certain circumstances, can be related to a spri ng-dashpot- 
hardware system. 
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2.2.4 Power Laws 

It can be shown that the use of discrete spring-dashpot models 
will allow us to obtain an accurate representation of viscoelastic 
phenomena over a large time scale provided we use a sufficient number 
of elements -typical ly one per decade. 

The idea behind the power law representation is to sacrifice 
accuracy at any one time in order to obtain a reasonable representation 
over the complete time scale between glassy (short term) behavior and 
rubbery (long time) behavior. This is often called broad-band approxi- 
mation. Landel and Fedors [38] gave a form for spectral repre- 
sentation H(x) guided by their understanding of polymer mechanics 

which may be written as, 

r \ n 


H(x) = C 



T 


exp (- t 0 /t) 


(2.14) 


v ) 


Clauser and Knauss [39] gave examples of this approximation of H(x) 
which are shown in Fig. 2.4. These results lead to the following 
broad-band approximation for the relaxation modulus: 


E r (t) = E 


+ 

e 



(2.15) 


The representation is seen to possess four arbitrary constants, 
including the limiting moduli values. The exponent n gives the slope 
of the relaxation curve (Fig. 2.5) through the transition region 
between glassy and rubbery behavior and x , for isothermal conditions, 
fixes the characteristic relaxation time. 
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Fig. 2.5 Relaxation Modulus E(t) vs Log Time 
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Recognizing that relaxation and retardation are basically 
reciprocal, one could also postulate a modified power law for the 
compliance spectrum, L(t), as 


L(t) = C' 



exp(-x/T 0 ) 


(2.16) 


However, the integration for the creep compliance does not give a 
simple form and it is more direct to assume a reasonable expression as 


D c (t) = D g + 



(2.17) 


where n may be the same exoonent as in the relaxation modulus and x' 

o 

is chosen for the best fit of the data. A special case results when 
t « x q and C = D e - D , 

n 


D c (t) = D g + C 



(2.18) 


A slightly modified form can be written as 

, n 


D c (t) * D g + C* 


(2.19) 


where x^ = A U, C* = C/A n , A is a dimensionless multiplication factor, 
and U = 1 is a time unit selected to accommodate data reduction (i.e. 

1 sec, 1 min, 1 hr, etc.). The power law in this form (Eq. 2.19) is 
often called the "Findley power law" and has been found to give an 
accurate representation of creep behavior of an extremely broad range 
of amorphous, crystal ine, crossl inked glass and carbon-fiber reinforced 
plastics and laminates. The equation has been verified mostly in 
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tensile tests, but in a few cases it has successfully described 
behavior in compression and under combined tensile and shear stresses 
as well . 

Upon conducting a screening experiment on 934 neat resin (at room 
temperature), we identified the exponent, n, in Eq. 2.19 as the most 
sensitive indicator of experimental errors. In addition, n was shown 
to be dependent on the duration of the creep test as shown in Fig. 2.6. 
Thus any short time creep test will result in an underestimation of n 
or in other words a stable evaluation of n is impossible as long as the 
flow term (C*t n ) is too small relative to the initial compliance D . 

The problem is easily handled, by performing a small additional experi- 
ment, which consists of recording the recovery strain after unloading. 
Creep recovery is mathematically given by 

€ r ■ 4«[(1 + X)" - ( A )"] (2.20) 

t — 1 1 

where A = — ^ — , as = a*C*t-| , and t-j is the time at unloading. 

Eq. 2.20 can be very easily derived using the power law for creep to- 
gether with the superposition principle. Fig. 2.7 shows that the n 
value obtained from recovery data is nearly insensitive to the duration 
of the recovery experiment. 

The creep experiment was repeated on T300/934 carbon epoxy with 
the load perpendicular to the fiber direction. The result plotted in 
Fig. 2.6 is compared with the result obtained on the neat 934 resin 
(Fig. 2.8) and shows that the composite n-value comes within 5% of 
the resin value after 3000 min. 
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Fig. 2.7 The variation of Power Law Exponent with 
Length of Recovery Test. 
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Fig. 2.8 Comparison of Power Law Exponent 
Variation with Length of Creep Test. 
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An alternative power-law type approach was developed by the 
Russian school, mainly under the influence of the ideas published by 
Yu. N. Rabotnov [23]. His formulation was based on the viscoelastic 
deformation of polymeric materials which can be described by an 
integral equation of the Bol tzmann-Vol terra type 


e(t) 


= £ [a(t) 


K(t - s)o(s)ds] 


( 2 . 21 ) 


where a(t) and e(t) are the stresses and strains in a uniaxial stress 
field, E is the instantaneous Young's modulus, and K(t) is the transient 
compliance or kernel function. Young's modulus can be determined by 
the basic quasi static measurements, while the kernel must include the 
smallest possible number of parameters that need to be determined 
experimentally. To satisfy the last condition, Rabotnov approximated 
the strain, e(t), with a power law and proceeded to deduce the kernel 
function K(t) which can be written as, 

°° Q n . n(l+a) 

K(t) = U “ J 0 rt(n+l)(l+a)J (2 ' 22) 

In (2.22), X, e and a are material parameters, r[(n+l )(l+a)] is the 
gamma function of the argument (n+l)(l+a), where 1 + a > 0. The kernel 
parameters can be determined graphically (when a full deformation curve 
is obtained) as well as numerically [19]. 

Functions of the type 2.22 are classified as exponential functions 
of fractional order. They behave as an Abel -type singularity for 
instantaneous loading and like an exponential for the succeeding period 
of time. A fairly detailed algebra for these operators has been 
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developed by Rabotnov and his collaborators. 

Equation (2.21) can be modified for the case of creep where a(t) 
is constant and becomes, 




K(s)ds] 


(2.23) 


For small t, the first term in (2.22) dominates, i.e., 

K < t >'HTCT (2 ' 24) 

By substituting (2.24) into (2.23), we obtain: 

6 <‘> 1 r n ♦ frkr tl+ “ ] < 2 - 25 > 

This again reduces to the Findley equation (2.19) 

D(t) = D Q + C*t n 

where C* = p (T^'n)' anc * n = I + a - It is clear again that Eq. 2.19 is 
only able to describe a limited portion of the full deformation curve. 
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Chapter 3 

SINGLE INTEGRAL REPRESENTATIONS BASED 
ON IRREVERSIBLE THERMODYNAMICS 


3.1 Introduction 

A single integral nonlinear viscoelastic constitutive equation 
has been developed by Schapery in a series of publications [34,40-44]. 
This chapter reviews his thermodynamic development with a special view 
to developing a better understanding of the nature of the four non- 
linearizing parameters which appear as well as their physical origin. 

The name irreversible thermodynamics is well established but 
confusing, since it suggests geometric irreversibility. This is not 
generally true because thermodynamic irreversibility relates to entropy 
production while a deformation accompanied by a change in entropy can 
be geometrically reversible as for example in rubber elasticity. Thus 
a name like "process thermodynamics" or non-equilibrium thermodynamics 
would be more related to the physics of a problem like viscoelastic 
deformation. 

The consistency of constitutive equations with the fundamental 
theorems of thermodynamics has been discussed by Truesdell and Toupin 
[45]. No restriction on constitutive equations is given by the first 
and the third fundamental theorems. The second excludes processes 
with negative entropy production and thus limits the field of constitu- 
tive equations, but not to the extent desired. It was demonstrated 
by Biot [46] that the more precise statements about the nature of the 
second law, which were made by Onsager [47], are sufficient to treat 
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problems in linear viscoelasticity. 


3.2 Gibbs Free Energy Formulation of the Schapery Equations 

This formulation is particularly useful, since it allows the 
treatment of temperature and stress tensor components as independent 
variables. The Gibbs free energy g is defined as 

g = u - Ts - - Q m q m (3.1 ) 

3 p M m tp 

where u is the specific internal energy, T the absolute temperature, 
s the specific entropy, p the density, q m (m = l,2,...,k) a set of k 
state variables consisting of observed variables (strains), and 
(m = l,2,...,k) is a set of k generalized forces defined by the virtual 
work condition S w = <sq.j . Eq. (3.1) shows that free energy can be 
accumulated by increasing internal energy by decreasing entropy or 
through a potential energy loss of the external loading. 

A mathematical expression for the strength of the entropy source 
as a result of "irreversibility" has to be derived first. We assume 
that the classical concept of entropy is extendable into non- 
equilibrium situations. This leads to the following result which is 
called the fundamental equation of the system 

Sj = S| (u,q.. ) (i = 1 ,2 , . . . , n) (3.2) 

The index I, identifies the system under study. The state variables 
q i (i = k+l,...,n) represent internal degrees of freedom.* These will 
allow us to describe micro-structural rearrangements related to 


*The number n of independent thermodynamic variables depends on 
the nature of the system. 
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straightening and relative sliding of long chain molecules. These 
variables are thus eventually observable, but certainly not control- 
lable. For this reason Q i (i = k+l,...,n) = 0 since any force 
associated with an internal degree of freedom would provide us with a 
means to control this internal coordinate. The entropy-increment 
dSj created by a change in du and n variations of dq^. is measured by 
the total derivative of the fundamental equation, i.e. 


dSj = 


as. 


3U 


n 

du + l 
q. i=l 


as. 


aq, 


u,q 


, d s- 


(3.3) 


where qj means all coordinates are held constant except q.. The incre- 
ment du is controlled by the first law of thermodynamics through 


1 

du = dh + l p~ Q. dq. 

i = l 1 1 


(3.4) 


Thus according to (3.4) any increment in internal energy is produced 
by an infinitesimal amount of heat (dh) crossing the boundary and/or 
an infinitesimal amount of mechanical work p~^ Q^. dq.. done on the system. 

The incremental process is reversible, when dq^ = 0. Then, 
du = dh and according to the second law. 


dh = TdSj (3.5) 

By substituting (3.5) into (3.3), we obtain 


In other words, the temperature measures the sensitivity of the 
internal energy for entropy changes. Thus (3.3), written now under 
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the restrictions of a reversible process becomes 


n 

Tds T = du + T 7 
1 i=i 


3Si 


3qi 


u,q 


, dq i 


(3.7) 


Since the partial derivative in the second term is an operation on a 

R 

state function Sj, it defines a state function , 

3s, 


-1 n R _ 
- p Q i - 


I 


(3.8) 


i 


Thus 


-1 „ R 


du = TdSj + l p' 1 Q i dq^. 


(3.9) 


i=l 


Eq. (3.9) is a result known in classical thermodynamics as Gibbs equa- 
tion and describes the associated reversible process, that causes the 
same increment of internal energy du as the irreversible process 
described by Eq. (3.4). 

By elimination of du between (3.9) and (3.4), we obtain 


Tds T = dh + l p" 1 (Q. - Q- R ) dq, (3.10) 

with Q. (k+1 , . . . ,n) = 0. 

Thus the entropy of a mass element changes with time for two 
reasons. First, because entropy flows into the volume element occupied 
by this mass element, second because there is an entropy source due 
to irreversible phenomena inside the volume element. This entropy 
source is always a non-negative quantity. 

Since the computational procedure needed for the derivation works 
only on isolated systems, we now create such a system by the artifice 
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of adjoining a large heat reservoir II to the physical system I. 

Eq. (3.11) states that the entropy is an extensive property, i.e. the 
entropy of the system is the sum of the entropies of the subsystems 


ds MI 

= dSj + dSjj 

(3.11) 

dS II = 

dh 
" T 

(3.12) 


Eq. (3.12) states that the reservoir plays the role of a reversible 
heat source, i.e. heat can be extracted from it without a change in its 
temperature T. A combination of (3.12), (3.11), and (3.10) leads to 


'I+II 


l n 

■sM, 


x i 


(3.13) 


D 

with X.j = Q-j - Qj • It turns out that the strength of the entropy 
source, Sj + jj has a very simple appearance. It is a sum of n terms, 
each being a product of a flux characterizing an irreversible process 
and a quantity, called thermodynamic force. 

In order to proceed from here we need to calculate the rate of 
approach to equilibrium. This is only possible by arguments outside 
the scope of continuum mechanics. The entropy source is explicitly re- 
lated to the various irreversible processes that occur in the system 
through a linear phenomenological relation between forces and fluxes 
q.j called Onsager's theorem. 


(. = y b. .(Q ,T)o. 

i >i ij'mT -j 


(3.14) 


This fundamental theorem states that when the flow q., corresponding 
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to the irreversible process j, is influenced by the force , cor- 
responding to the irreversible process i, then flow is then also 
influenced by the force X. through the same influence coefficients, 

J 

b..(Q ,T). Eq. (3.14) can be rewritten as 
ij m 


Q. = Q. R + b. .(Q ,T)q . 


(3.15) 


It is shown in Appendix B that: 


n rc _ 

Q,- = P 


9S_ 

aq. 


i = k+1 ,. . . ,n 


The generalized force Q. associated with q^ through the Gibbs equa- 
tion (3.9), is derivable from the Gibbs free energy g (Q m , q r , T) 


n m 


- p 


aQ m 

m 


m = 1 ,2 , . . . ,k 


The problem is thus reduced to the solution of two systems of equations, 


q m 


32 _ 

aQ m 

m 


m 1 ,2, . . . ,k 


(3.16a) 


p I s - + b (Q ,T)q = 0 
aq r rs vx m’ /M s 


r,s = k+1 , . . . ,n 


(3.16b) 


A constitutive assumption on the dependence of the Gibbs free energy 
on Q m , q^ and T has to be made. This is accomplished by a Taylor series 
expansion around the equilibrium state such that, 

pq = c n + c 0 + c q + i d Q Q + d Q q 
y R m 4 m r H r 2 mn v m w n rm y m t 


+ 7 d rs q r q < 


(3.17) 


where c^, c- and d. . may depend on temperature, but not on q r and Q m< 
Note that there exists a distinct relation between summation index 
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and physical dimension for these coefficients, i.e. 


Coefficient 

Represents 

Summation Range 

C R 

Energy 


c m 

Mechanical strain 

m = 1 ,2, . . . ,k 

c r 

Thermodynamic force 

r = k+1 , . . . ,n 

^mn 

Mechanical compliance 

n,m = 1 ,2,. . . ,k 

d n» 

Transforms hidden variables 
g r into mechanical strain 
units 

r = k+1 , . . . ,n ; 
m = 1 ,2,. . . ,k 

d rs 

Thermodynamic modulus 

r,s = k+1 ,. . . ,n 


A further approximation of the free energy is found by an expansion 
with respect to the uncontrollable variables only, i.e. we can be 
far away from equilibrium with respect to Q . 

09 = 09 R + c r q r + ? d rs q r q s (3J8) 

g D , c' and d' can be functions of the independent variables Q m and T. 
k r r s in 

Eq. (3.18) must reduce to (3.17) under sufficiently small forces, 
which leads to the relations 


5dn = Cr,+ C Q + 7T d Q Q 

y R R m v m 2 mn ^m 4 n 


(3.19) 

(3.20) 

(3.21) 

The various terms in equation (3.18) have the following physical 
significance; 


c' = c + d Q 
r r rm *m 


d' = d 
rs rs 
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• pg r An energy term that acts as a potential for the time inde- 
pendent response. The dependence of c and d on stress 

m mn 

and temperature can be lumped together with c^ into a 
single nonlinearizing parameter. 

•c'q = ( c + d Q)q = $ 0 q +d Q q This term is de- 
r M r v r rm w rrr H r r H r rm y m M r 

fined by the energy associated with a simultaneous action 

on the internal coordinates q^ and the controllable 

r 


variables CL and 0. 
m 


1 


• j d^, s q r q $ This term is defined by the interaction energy 
associated with the simultaneous action on internal co- 
ordinates q r and q . 

The above items which make up equation (3.18) are of particular interest 
as they define the basic interaction mechanisms involved in three dif- 
ferent nonlinearizing functions which appear in Schapery's final single 
integral equation as developed subsequently. 

Substitution of (3.18) into (3.16b) gives 

dq c 


c' + d 1 (Q ,T)q + b (Q ,T) -r*- 
r rs vy m /H s rs vy m’ ' dt 


(3.22) 


In order to reduce (3.22) to a system of n-k differential equations 
with constant coefficients, the further assumption is made that, 


b rs<V T > * W T > b 


rs 


(3.23) 


in which b r$ are now defined as viscosities. This equation implies that 
the flow is Newtonian but with stress and temperature dependent 
viscosity. The implicit stress and temperature dependence can be made 
explicit through a shift function a^Q^T)- A 017 experimental 
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evidence exists for the special case where all viscosities have a 
common shift factor independent of stress. In such cases a Q becomes 
the familiar a-j-(T) shift factor found in polymer literature [28]. 

The following form for dj, is guided by (3.23), 


r,s = k+1 , . . . ,n 

jR 


(3.24) 


d; s (Q m> T ) - a G (Q m ,T)d" s 

where the constant symmetric "thermodynamic" modulus, d^ $ , is now de- 
fined at the reference temperature T R and Q m = 0 as indicated schematical- 
ly in Fig. 3.1 and a g(Q m >"0 is a shift function as defined below. 

< * d ;s<°' T R> ’ d rs< T R> < 3 ' 25 > 

It follows from Eq. (3.13), (3.14) and (3.23) that ag nonl inearizes 
the entropy production. 


S I+II ll a o( Q m’ T) b ij ^i 


(3.26) 


On the other hand a^ nonl inearizes the third term in the Gibbs free 
energy expansion (3.18) and therefore pg is proportional to the above 
quantities as indicated by, 


pg - [ a s (Q m> T ) d R rs q r q s (3.27) 

We cannot make an a-priori statement on the dependence of a^ with 
stress, but the dependence of d p S (Q m ) at a certain temperature T can 
increase, decrease, or be independent of the force Q m (Fig. 3.2). 

A constitutive assumption on a^ has been developed by Schapery 
[48] for the case where microcracking and flaw growth (in the opening 
mode) are the dominant softening (damage) mechanisms. The reciprocal, 
l/a^, is an increasing function of the q-th order Lebesque norm of 
tensile stress, <t>', as specified by, 
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(3.28) 


±~ = 1 + 




G(4>) d$ 


where , 


l _ 


a q d£ 


(3.29) 


G(<t>) is positive and reflects the statistical distribution of damage. 

r 

The quantity, G( 4 > ) d4> , can be found from mechanical property tests 

V 

wherein q is a positive constant. 

Extensions of Eq. (3.29) have been made for multiaxial stress 
conditions through incorporation of the first and second invariants of 
the stress tensor (e and J 2 ), while aging and/or rehealing is intro- 
duced through a function M(t) such that. 


I _ 




M n [A i9 + A 2 /J^ ] q dC q 


(3.30) 


Here, A^ and A 2 are constants which satisfy the relation A 2 = (1 - 

A-j ) and when they are both positive there is a suppression of crack 
growth by (9 < 0). 

Due to the simplifications introduced in (3.23) and (3.24), a 
system of n-k differential equations with constant coefficients 
solvable for q^-j , . . . ,q n now remains. These together with the k 
algebraic equations obtained by substituting the free energy expression 
(Eq. 3.18) into the k equations of (3.16a) are. 


a G d rs "s + a D b rs 3T = ’ ^ 


dg R 


q m = " p 30 


m 


K a 

Q q r 
m 


(n-k equations) (3.31) 

(k-equations ) (3.32) 
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The former equations (3.32) can be rewritten with the reduced 


time tjj as. 


d 


R 

rs 


dq s 

+ b rs3T 



(3.33) 


where, 

d4> = dt/a o and a a = a o I a 6 

The final system of equations can be explicitly written as, 


3c k+l 

3Q, 

. 8c k+2 
3 O'! 

3C n 
* 3Q 1 

0 

0 

0 


+ 

cr 


q l 


3g R 

3C k+l 

3Q 2 

3c k+2 

3Q 2 

K 

" 3Q 2 

0 

0 

0 


q k+2 


q 2 


0 aq^ 

3c k+l 

3C k+2 

3Q k 

K 

■ - 3Q k 

0 

0 

0 


% 


q k 


ag R 

h r 

Vl ,k+l 

h r 

°k+1 ,k+2 ' 

d R 

d k+1 ,n 

b k+l ,k+1 

b k+1,k+2 

• • b k+l ,n 


d Vi 

cty 


• C k + 1 


0 

h r 

V2,kt1 

h r 

k+ 2 ,k+2 

d k+2,n 

b k+2,k+l 

b k+2 ,k+2 • 

• b k+2 ,n 


dq k+2 

dii» 


' c k+2 


0 

d n,k+l 

h r 

n,k+2 

.. d R 

n, n 

b n,k+l 

b n,k+2 

b n,n 


dq n 

w 


- c n 


0 


(3.34) 

Remark that the matrix partition on the upper left hand corner side 
has a maximum dimension of 6 x (n-6) while the two lower partitions 
are square matrices with dimensions (n-6) x (n-6) where n can be very 
large, since this is the number of internal coordinates, needed for an 
adequate description of the inelastic material phenomena. 
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p 

It can be proven that a simultaneous di agonal ization of d£ s and 
b rs is always possible, which means physically that a transformed 
space exists where the n-k internal deformation modes do not interact. 
Thus, 


b 


rs 



and 



for 


r = s 
r f s 


(3.35) 


and 


dq r 

d q + b -r^ = - - 
r M r r dip c 

The solution of (3.36) is 


c ; 


q r - (1 - e~* hr ) 


f 

1 

fi) 

w 

. a G 


(3.36) 


(3.37) 


with x = b r /d r . Eq. (3.37) has a fading-memory Kelvin-type kernel and 
can thus be seen as the evolution equation for a Kelvin unit associated 
with the physical process represented by q^. A point of interest here 
is that there is no force Q r associated with the (n-k) functions q r 
which describe the evolution of the internal system. It appears however 
that these internal mechanisms are driven by forces - c^/a^ spiked 
into the Kelvin unit through a force transfer mechanism, which was 
mathematically given as and includes a thermodynamic force c^ due to 
temperature differences (<t> = T - T R ) . 


c' = c + d (Q ,T)Q 
r r mr y m’ /y m 


(3.20) 


The transfer function (Eq. 3.20) describes how forces Q m associated 
with the observed coordinates (q-j ,q 2 « . . . .q^) diffuse down into the 
internal material structure where they drive the hidden coordinates 
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This force transfer process is not necessarily linear and introduces 


a third nonlinear function Q m through the equation 

c ' =3 d> + d R Q (Q ,T) 
r p r T rm w nr y n’ ‘ 


(3.38) 


where ct> is T - T D and s is the thermal modulus. 
k r 

Note that functions d^ s (which show up in the third term of the 
free energy expansion (3.18)) and d^ (which show up in the second term) 
represent two completely different nonlinear-physical phenomena. This 
reasoning leads to the two distinct functions a^ and Q m respectively. 

Substitution of (3.37) into (3.31) gives 


'’m 


9g 

p 3Q 


l 


3C 


m 


sQ 

r m 


(1 - e 


-tyhr) 1 


d r a G 


(3.39) 


and for time varying forces. 


g m = 


9g R . r 9c r 1 
P ^m r 3 ^m d r J 


(1 - e 


-(♦-♦•>A r) d(c r /a G > 


w 


d<p‘ (3.40) 


By substituting (3.38) in (3.40), 
3g R , Sq n 


q m = - P 
m 




d(Q n /a r ) 

AS (ip-ip ‘ ) dx 


n P 


Acx^ ( ip— ip ' ) 


dU/a r ) 

' 


dx 


dx 


(3.41) 


where 


AS = l (1 

np £ np 

(r) 


e ^ /x r) 


Aa n - l 0 - 


(3.42a) 

(3.42b) 
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(3.42c) 



sQ n 

9 1 = aQ7 


(3.45b) 

(3.45c) 



(3.45d) 


By way of summary, the nonlinear parameters g , g-j , g 2 and a Q 

have the following physical significance. 

• q Nonlinearizes the elastic response and can in principle 
o 

increase or decrease with stress. The former is described 
as a softening phenomenon while the latter is considered 
to be a hardening phenomenon. 
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• g-| Experimental results indicate that g-j is an increasing 

function of stress which requires Q-| to be a monotonically 
increasing function of Q-| as seen in Fig. 3.3. This means 
that the stresses that are transferred to the microstructure 
may have a higher intensity than the external force. Note 
also that when the material is linear viscoelastic, 
g-| = 1 and Q-j is a linear function of Q-| . 



Fig. 3.3 Dependence of Generalized Force Q-| or Q-j 


• g£ This factor contains two nonlinearizing functions, Q-j and 
a^. It can be seen as the ratio of the force that could 
be spiked into Kelvin-type microunits taking into account 


that the moduli of the springs in these units change by a 
factor a^. This change can be stiffening (ag > 1) or 
softeni ng (a^ < 1 ) . 


* at 

a 



Represents a stress and temperature dependent shift 
along the time axis. Here a D (Q m ,T) decreases with 
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temperature and with increases in stress-induced free 


volume. Also, a g(Q m »T) decreases with temperature 
but an increase in stress would produce more shift. 

The Kelvin type kernel aD in Eq. (3.44) can be approximated by a 
power law which was proved to be valid under restricted conditions in 
Chapter 2, 

aDU-*') = l D (1 - e'^ ,)/T n (3.46) 

r 

or 

A D ( ip-tp ' ) ~ c ( ip-tp ' ) n (3.47) 

With this approximation Eq. (3.44) becomes. 


e = g 0 D q a + g. 


t n d(g ? o) 

C ( tp-ip ' ) — — dx 

o 


(3.48) 


For a creep test a = a Q H(t), eq. (3.48) reduces to, 


6 = 9 0 D o °0 + 9 1 «2 C 




(3.49) 


The experimental procedure that should be followed to obtain the 
three linear parameters (D ,C,n) and the four nonlinearizing 
parameters ( g Q ,g-j > g 2 » a CT ) has been described in detail by Lou and 
Schapery [49] on the basis of a graphical shifting procedure based on 
creep and creep-recovery data. A computerized procedure that avoids 
tedious graphical shifting was developed during the current research 
effort by Bertolotti et al . [50]. 

The general equation (3.44) can be simplified to model specific 
thermorheological complex materials. For sufficiently low stresses, 

- Q m = Q m = a which means that g-j = 1 and g 2 = Va G - Also, a Q and a D 
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are only functions of temperature (ap/ag = a-p). Eq. (3.44) now reduces 


to 


rt 


e = g o a + 
J o o 


AD(ip-ip ' ) 


( \ 
a 




dx 


(3.50) 


where log ^ = log t - log a-p. For a = a Q H(t) this becomes 


D t = Dj(T) + AD(ip)/a G (T) (3.51) 

with Dp = e | o Q and Dj(T) = g Q (T) D Q . Thus, 

log (Dj - Dj) = log aD - log a^ (3.52) 


Eq. (3.52) shows that master curves can be obtained by rigid horizontal 
| log a-p| and vertical 1 1 og a^ | shifts. 

Eq. (3.52) suggest that the machinery needed for accelerated 
characterization can only be effective when the vertical shift ag(T) 
can be separated from the horizontal shift a-p(T). This is only 
possible with additional transient temperature tests as suggested 
from Eq. (3.41). Some actual data were obtained by Schapery and Martin 
[51] and very recently by Weitsman [52]. In the former, a theoretical 
form for ag(T) which applies below and in the neighborhood of the 
glass transition temperature Tg was given as. 


a = y in (1 + z) 

6 exp [ y 2 /(l + z)] - 1 


(3.53) 


where y = aH/RT r and z = (T - T R )/T R . Eq. (3.53) is based on kinetic 
theory wherein aH is the constant activation enthalpy for hole forma- 
tion. Simple forms for ag on the basis of rubber elasticity are well 
established at and above the Tg. 
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The 

materials 


results discussed herein have been applied to composite 
by Schapery and his co-workers [49,53]. 
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Chapter 4 
DELAYED FAILURE 


4.1 Introduction 

Failure is most often treated as a separate issue from the deter- 
mination of modulus properties of materials. In fact, most failure laws 
are derived empirically from observations related to a catastrophic 
event such as yielding or rupture. As a result, a great deal of 
testing and data analysis is necessary to establish an appropriate 
failure law for a material. On the other hand, modulus or constitutive 
laws are derived by more rational means of relating deformations to the 
forces which produce them. For this reason, often much less testing is 
necessary to define a constitutive law for a material especially if 
deformations do not depart from the elastic or reversable deformation 
range for a material . 

Failure, however it is defined, should be a part of a complete 
constitutive description of a material. In other words, the key to 
dealing effectively with failure lies in treating its behavior as a 
termination of a nonlinear viscoelastic process. Perhaps, for this 
reason, a number of investigators have suggested that modulus and 
strength laws should be related to each other for viscoelastic materials 
[4,38,54]. 

The concept of distortional energy as a measure for the criticality 
of a given state of stress in an elastic material (von Mises criterion) 
cannot be carried over directly as a governing criterion for failure 
in viscoelastic materials, mainly because viscoelastic deformation 
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involves dissipative mechanisms. Thus at any time the energy balance 
can be written as: 

Total deformation = Free (stored) + Dissipated 
energy energy energy 

In case of a perfectly elastic/viscoplastic material, inelastic 
deformations occur only after a limiting value of stress has been 
exceeded. Owing to this assumption, such a yield criterion does not 
differ from the known criteria of the inviscid plasticity theory. The 
most appropriate form on the latter is often based on the energy approach 
(Huber-Henkey-von Mises) in which a critical value of the conserved 
(accumulated) distortion energy defines the transition from the elastic 
to the inelastic state. 

In case of an elastic-viscoplastic material, the rheological 
response occurs from the beginning of the loading process. Under such 
conditions it seems appropriate to take into account the accumulated 
energy W^, as well as the dissipated power where both quantities are 
regarded to be responsible for the behavior of the viscoelastic material. 
The yield condition is then given by a function with two arguments 
such that 

f(W E ,W D ) = o 

where the specific form is based on experimental evidence. The simplest 
possible form is the linear combination 

5 E W E + 5 D C ^D = k 

with 

0 < U E >C D } 1 1 
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With 5 D = 0 we assume that the dissipated energy has no effect on 
the yield criterion. This corresponds to Reiner and Weissenberg 's 
point of view. Reiner and Weissenberg [55] suggest that the energy 
storage capacity is responsible for the transition from viscoelastic 
response to yield in ductile materials or fracture in brittle ones. 

They assume that a threshold value of the distortional free energy, 
called the resilience of the material, is the governing quantity. It 
is clear that the Reiner-Weissenberg criterion applied to a material 
with zero dissipation (elastic material) becomes identical to the von 
Mises criterion. When applied to a viscoelastic material, however, we 
have to keep track of the free energy which is history dependent. If 
the mechanisms through which total deformation energy is transformed 
into dissipated energy are so activated that no free energy can be 
accumulated, there is practically no limit to the amount of deformation 
energy which can be applied without impairing the strength. That is, 
forces up to a certain magnitude can be applied for any length of time 
without leading to rupture. If it turns out that the material cannot 
accommodate this energy redistribution fast enough then the material 
will store energy initially but after a period of time failure will 
occur. The failure process is therefore delayed which limits the life 
of a given structure to a finite value. The instant of yielding is 
thus clearly dependent on the final outcome of a conjection between 
deviatoric free energy and dissipated energy. Thus the effect of the 
strain history on the delayed yielding phenomenon follows from this 
model in a natural way. In this respect, it should be noted that both 
Naghdi and Murch [56] and Crochet [57] have used a modified von Mises 
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criteria for viscoelastic materials by assuming, ab initio, that the 
radius of the yield surface Y does depend upon the strain history, 
through a Euclidian metric, x> in strain space. Their formulation 
can be expressed as, 


f ' 7 hj s ij - y2 * 0 • ¥ * Y <*> 


x = [(«Y. _ e !.)(eY. - e ^.)] 1/2 

* LV ij ij /v ij ij' J 


(4.1) 

(4.2) 


where x represents the distance between the purely elastic strain 

response e*r. and the strain state corresponding to yield eY. as is 
' J *1 J 

geometrically shown in Fig. 4.1. The functional relationship between 



Fig. 4.1 Geometric interpretation of the Euclidian metric x (as de- 
fined by Crochet). 

Y and x for one-dimensional behavior was given by Crochet [57] to be, 

Y(t) = A + Be' Cx (4.3) 

where A, B and C are material parameters. The radius of the yield sur- 
face is thus a continuously decreasing function of the "distance" between 
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two states of strain (e E . and e v .). 

* J ^ J 

An extension of the Naghdi and Murch-Crochet approach to fiber 

reinforced material systems is possible when we can relate the stress 

in the matrix phase (a) to the applied external stress {a}. The 

m 

introduction of stress concentration factors, which relate both these 
stress-fields is suggested by Hill [58] 


{a} = B {a} 

m m 

and equation (4.1) becomes 


(4.4) 


f = {o} T [B m ]{a} - Y 2 = 0 (4.5) 

in which Y is the tension yield stress in the matrix, as defined by 
eq. (4.3). It should be noted, however, that eq. (4.5) only describes 
initial yielding. The incorporation of a hardening rule, together with 
a flow rule should be considered in order to model accurately the 
phenomena after initial yielding up to final fracture. 

The extension of the Reiner-Weissenberg approach to a general 
anisotropic material system can be done quite easily when we use the 
general definition of distortional energy, as defined by von Mises [59] 


W dist 


2 ^ S ll a x + s 22 °y + + S 66 T z^ + s 12 a x CT y + 


+ S 1 4 Vx + + s 66 T y T z 

1 ^x^ll + S 1 2 + S 1 3^ + ••• t z^ s 61 + s 62 + 


63' 


’ll 


+ S 22 + s 33 + 2 ( s i ^ + s„ + s-,-, ) 


12 3 23 


31 


(4.6) 


The advantage of this approach is that the onset of yielding or failure 
is controlled by only one material constant st , while the former 
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approach (Naghdi et al .<•) introduces at least three parameters. It has 
been argued by Pagano [60], on the other hand, that the distortional 
energy, as defined by von Mises, has only physical significance in the 
special case of an isotropic material. The full implications of 
Pagano 's comments are unknown at the present time. A possible alter- 
native way to extend the Reiner-Weissenberg criterion, would be to 
evaluate the free energy in the matrix phase. This could also be done 
through Hill's concentration factors. 

The following is a brief review of the Reiner-Wei ssenberg criterion 
which follows Bruller's [61-68] extensive investigation along these 
1 i nes . 

4.1.1 Free Energy Accumulation for a Linear Kelvin Chain in Creep 

For a single Kelvin unit in series with a spring as shown in 
Fig. 4.2, 



Fig. 4.2 Linear Kelvin model 
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€ 


(4.7) 


= € g + e l 


where 


€ 1 = D 1 ' exp (' t/x l ) ] 

or 

e(t) = [D g + D^l - exp (-t/x-j)]^ 

The stress, a , in the spring with compliance D-j is, 
e . 

= q^- = a Q [1 - exp ( ~t/ T-| )] 


(4.8) 


(4.9) 


(4.10) 


The stress, a^, in the dashpot is then, 

°d = CT o exp ^' t/x l^ (4-11) 

A physical interpretation of free energy would be the energy stored in 
the springs D g and D-j , 

? * D n ft 

W . = a + o 6 dt 

springs o 2 j s s 

o D _ ? D 1 O 

= a o 2^ + °o T [1 " exp (- t/x l^ (4.12) 

In the same way we can calculate the energy dissipation in the dashpot. 


W 


dashpot 


o d e d dt 


2 D 1 

“ a o r - ci " exp (' 2t/x i 


(4.13) 


The total energy in the three parameter system is given by. 
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w = w + w 

total springs dashpot 


= a. 


U 

^ + D 1 (1 - exp ( - 1 / t i ) 


(4.14) 


The extension of Eq. (4.12), (4.13) and (4.14) to a n-element Kelvin 
unit with a free spring is straightforward. 


W -a 

springs o 


D N D. 

** + l (1 - exp (-t/x.))‘ 
L 1=1 


N D, 


W dashpots = CT o T (1 ‘ exp ( - 2t/T i )} 


W total °o 


'D N 

p 3 - + l D. (1 - exp (-t/x- ) ) 
c i=l 


By taking the limits for t -* « of these equations, we obtain 


(4.15) 


(4.16) 


(4.17) 


Springs a o 


"D N D." 

_£L + V _L 

2 i=l 2 


2 N D i 
W dashpots = °o 2" 


(4.18) 


(4.19) 


J im W total = % 
t-*» 


D N 


(4.20) 


Thus half the work done by the external forces goes to increase the 
free energy while the other half is dissipated. 


4.1.2 Power Law Approximations for Free Energy 

The free energy expression in the form of Eq. (4.15) forces us to 
identify the 2N material constants D.. and . This can be avoided 
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using a power law approximation as follows, 


W 


total 


D N 

' ^ ' l 1 .! "■ ° 


exp ( - t/ t i ) ) - Dt n 


(4.21) 


and 


W 


dashpots 
2 
a o 


N 


j l D. (1 - exp ( — 2 1 / t - ) ) - i D(2t) n (4.22) 

i=l 1 


This allows us to obtain the following simple expression for the free 
energy 


W = a 2 

springs o 


'fail ure 


2^- + D (t n - j (2t) n ) 
/w 


(4.23) 

(4.24) 


j 2 + D (t n - \ (2t) n ) 

The nonlinearizing parameters introduced in Eq. (4.24) lead to 
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(4.25) 


The ability of Eqs. (4.15) and (4.25) to model experimental data is 
checked in Chapter 7. 


Note on the Zhurkov Criterion 

The existence of a strength limit as a real physical characteris- 
tic of a material was seriously doubted by Zhurkov [69] due to experi- 
mental evidence that solids may also fracture at stresses below the 
strength limit. He set up systematic and careful measurements of the 
lifetime of solids under a constant tensile stress. Both tensile 
stress and temperature were widely varied among experiments. More than 
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a hundred substances of 'all principal types of solids were investigated-- 
metals, alloys, glasses, polymers , crystals. It turned out that the 
time to fracture possessed a surprisingly uniform dependence when 
suitably expressed in terms of the stress (a) and temperature (T). 
Analytical treatment of this vast amount of data led him to the uni- 
versal formula for the temperature and stress dependence of lifetime, 
t, (except for very small stresses), 

t = t q exp [(U Q - ya)/kT] (4.26) 

where k is the Boltzmann constant, t q is a constant on the order of 

- 1 3 

the molecular oscillation period of 10 sec., II is a constant for 

o 

each substance regardless of its structure and treatment, and y depends 
on the previous treatment of the substance and varies over a wide 
range for different materials. 

Equation (4.26) can be formally incorporated into a continuum 
mechanics analysis in a modified form to account for time varying 
stresses and/or temperatures, 

•X dt 

T o exp L(U q - ya(t) )/kT (t) J " 1 (4.27) 

Roylance and Wang [70] used Eq. (4.27) in a finite element code where 
the current value of the above integral was computed at each node. The 
time and location of fracture was determined when the integral value 
reached unity at any node. 

Williams and his collaborators [71,72] disagree with Eq. (4.26). 
They argue that the rupture process is much more detailed than re- 
flected by this equation. Their experiments on stress history induced 
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rupture led them to the conclusion that "Zhurkov's impressive findings 
appear to be the somewhat fortuitous result of a particular test 
method. " 

When Eq. (4.27) is applied to such highly heterogeneous material 
systems as composites, the question of whether or net fracture can be 
regarded as a thermally activated process of damage accumulation is far 
from settled. It is also obvious that x , U Q and y lose their physical 
meaning. Regel et al . [73] suggest a life of mixture type predictions 
for these variables. 
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Chapter 5 

RESIN CHARACTERIZATION: EXPERIMENTAL RESULTS 

AND DISCUSSION 


5.1 Introduction 

In order to understand experimentally observed discolorations of 
the Fiberite 934 neat resin due to postcuring and some anomalities of 
the creep response in the temperature region, 200-250°F (93.3-121 °C) , 
it was decided to study details of the relationship between resin 
chemistry, various postcuring procedures, and mechanical properties. 
General resin information is summarized in 5.2 while experiments and 
analysis related to the mechanical characterization of the resin are 
reported in 5.3. 

5.2 General Resin Information 

5.2.1 Resin Chemistry and Curing Procedure 

Many composites, especially carbon-fiber based composites, utilize 
an epoxy resin matrix. From the hundreds of different epoxies that 
are available, only a small number (~ 10) are used for composites. The 
epoxy of interest in this study is a polyfunctional epoxide based on 
the tetraglycidyl derivative of methyl enedi aniline (TGDDM), available 
from Ciba-Geigy under the tradename MY 720. This epoxide, which is now 
the basis of virtually all aerospace matrices for carbon fibre both 
here and abroad, is shown in Fig. 5.1a. The Ciba-Geigy Eporal hardener 
used in conjunction with MY 720 is 4,4' diaminodiphenyl sulfane (DDS) 
Ciba-Geigy Eporal and is shown in Fig. 5.1b. 
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a) Tetraglycidyl 4,4' di ami nodi phenyl methane epoxy TGDDM 
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o 

b) 4,4' di ami nodi phenyl sulfane (DOS) 

Fig. 5.1 Chemical Structure for Epoxy Resins 

Because the DDS cure is sluggish, even at moderately high tempera- 
tures, a boron tri fluoride catalyst is used ( B ) in order to produce 
a more manageable cure cycle. Despite the fact that these ionic 
curing agents often have poor humidity resistance, they are used as a 
compromise between mechanical properties and manufacturing convenience. 
The final epoxy network structure depends not only on the chemistry of 
the system, but also on the initial epoxy versus curing agent ratio 
and cure conditions. At the beginning of this research program, we ob- 
tained three different panels manufactured with three different curing 
procedures. Cure details are given in Appendix A. 
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It is usually assumed that the properties of the cured resin in 
a composite are the same as those of the bulk material. However, 
McGullough [74], Hancox [75] and Yang, Carlsson and Sternstein [76] 
suggest otherwise. The reasons mentioned by these authors are: 


• The structure in a small volume of resin may be more ordered 
than in the bulk material, because the exotherm in the 
composite will be limited to the overall curing temperature 
by the presence of good conducting fibers. 

• The small distance between the fibers may provide steric 
hindrance and inhibit epoxy-amine cross-link reactions 
during the latter stages of cure and limit the overall 
achievable cross-link density. 

• Due to the large amount of surface area of the fibers, 
adsorption of epoxy segments at the fiber surface would 
cause a tighter network. 

• A state of negative hydrostatic stress could be imposed on 
the epoxy due to the fibers, restricting the volume expan- 
sion of epoxy during heating. 


No convincing evidence has been produced yet (to this author's knowledge) 
that one or a combination of these mechanisms leads to the resin in the 
composite that is thermorheologically different from the neat resin. 

As a polymer is cooled at a finite rate through the glass transi- 
tion into the glassy state, a sudden increase in modulus parallels a 
rapid decrease in molecular mobility. The process of molecular re- 
arrangements operates on its own internal time scale. Equilibrium can 
only be achieved if the cooling rate is so slow that the molecular re- 
arrangements are given sufficient time to operate. If the cooling rate 
is too fast, these rearrangements are essentially quenched into a non- 
equilibrium state, i.e., a state characterized by excess entropy, 
enthalpy and volume. This gives rise to a driving force towards 
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equilibrium, i.e., the system densities till the free volume, excess 
entropy and excess enthalpy all reach zero. This process can operate 
on a scale of days, months or years and is known as aging. 

Experimental equilibrium has been unexpectedly difficult to 
achieve. Spontaneous changes in the properties of glassy polymers have 
been observed long after vitrification. 

Excess entropy cannot be measured directly but excess enthalpy 
is easily measured by calorimetric methods (DSC). The latter show a 
decrease in enthalpy with time in a similar fashion as the time dependent 
mechanical properties as reported by Petrie [77,78], As indicated in 
Fig. 5.2, she did not find a one to one correspondence to specific 
volume behavior. 

In our experiments, we tried to achieve this equilibrium state by 
a postcuring of 350°F ± 10°F (176°C±5°C) for 4 hr ± 15 min. followed 
by a slow, controlled cooldown at a rate of 5°F/hr (2 .5°C/hr)(see Fig. 
5.3). 

5.2.2 Thermal Transitions in the Epoxy System 

Two distinct transition temperatures at 140°F (60°C) and 356°F 
(180°C) respectively have been identified by Yeow and Brinson [4] dur- 
ing their experiments on T300/934 graphite epoxy. The latter was the 
primary glass-transition temperature, T^, while the former was a 
secondary transition. Their measurements were based on detecting dis- 
continuities in the coefficient of thermal expansion which were measured 
using electrical strain gages. 
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Temperature ( deg F ) 


o 



Annealing Time ( hr) 

Fig. 5.2 Shear Storage Modulus, G-j , Logarithmic Decrement, a, and 
Enthalpy Relaxation, AH, of Quenched Atactic PS vs. 
Annealing Time at 92°C. G-J and A Were Measured at 0.6 Hz. 
Petrie [77,78]. 1 



Time ( hr ) 


Fig. 5.3 Postcure Cycle Used for 934-Resin and for T300/934 Composite. 
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Keenan, Seferis and Quilivan [79] used oscillating mechanical 
fields to characterize the dynamic mechanical behavior of an epoxy 
system. Their results are shown in Fig. 5.4 with an aj-transi tion 
centered around 212°F (100°C) and an a-transition at 464°F (240°C). 
However, these definitions of Tg suffer from the fact that they are 
based on dynamic mechanical quantities which are quite sensitive to the 
frequency of experimentation as well as the sample heating rate.* Con- 
sequently, the Tg values defined from dynamic mechanical experiments 
are invariably higher than those obtained by a more conventional 
method, i.e., dilatomeric experiments. 



Temperature (°C) 

a) Storage Modulus 



Temperature (°C) 

b) Loss Modulus 


Fig. 5.4 Storage Modulus and Loss Modulus after Keenan et al . [79]. 


*The conclusions obtained by Keenan et al . can be carried over 
to the 934 resin, since this particular member of the TGDDM-DDS 
epoxy family has the same basic network. Any difference should be 
one of degree, not in kind. 
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5.2.3 Environmental Stability 

1 ) Thermal Stability 

Different IR analysis of an epoxy postcured above 350°F (180°C) 
reveals an additional broad sorption band in the 5.5 - 14.3 y region 
according to Morgan [80]. The intensity of this sorption increased with 
exposure time and increasing temperature. A strong IR sorption at 
5.85 y appears during the initial stages of the degradation process, 
which indicates the formation of carbonyl groups. Such groups should 
result from oxidation of unreacted epoxide rings to form a-hydroxy 
alolehyde and carboxylic acid [81]. 

Degradation causes a discoloration of the 934 epoxy as indicated 
in Fig. 5.5 and may indicate a significant change of behavior in the 
glassy region due to postcuring. We observed that heating above 356°F 
(180°C) for 36 hr seemed to cause embrittlement, i.e. decreased tensile 
strength and ultimate elongation accompanied by a slightly increased 
modulus (Fig. 5.6). In addition, an increased moisture absorption 
rate (Fig. 5.7) was found. This is most likely due to a thermal oxi- 
dation phenomenon. Figure 5.8 shows the storage modulus E' and the 
damping ratio 6 for the same resin. 

2) Moisture Absorption 

One consequence of diffused water is that it lowers the mechanical 
rigidity of the resin. In extreme cases the effect is to lower the 
glass transition temperature by as much as 212°F (100°C). The effect 
of moisture on the glass transition temperature and the rates of 
absorption and desorption for different resins (including the Fiberite 
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4 hours at 350° F ( I77°C), 
slow cooled , 5° F 
( 2.5 ° C ) per hour . 


15 minutes at approximately 
10-15° F ( 5- 8° C ) above T g , 
then quenched down to room 
temperature. 


Short overheating at 500 °F 
( 260 ° C ) , then quenched down 
to room temperature. 

| 


36 hours at 380° F (I93°C), 
slow cooled , 5 ° F 
(2.5° C) per hour. 






Fig 5.5 Effect of Cure Cycle on Resin Color 
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5.7 Effect of Postcuring on Moisture Absorption Rate of 934 
Resin at Room Temperature. 














934) has been studied by Deisai and Whiteside [82]. Figure 5.9 shows 
the relationship between equilibrium moisture concentration and primary 
glass transition temperature for several epoxy resins. Figure 5.10 
shows the effect of the humidity in the test environment on the 
equilibrium moisture content in the same resins. Figure 5.10 shows 
that a prolonged postcure period could adversely affect the moisture 
absorption rates, indicating a damaged matrix. The specimens were 
allowed to pick up moisture in the lab [75-80°F (24-27°C) and ~ 45% RH] . 
Figure 5.11 shows the influence of moisture on the w-transi tion . It 
is seen that moisture has a strong effect especially on the width of the 
w-transition. 



Moisture Content, %(wt) 

Fig. 5.9 Effect of Moisture Content on [82], 
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Fig. 5.10 Effect of Relative Humidii 
tration in Epoxy Resin 
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The specimens we used for creep and creep-recovery lost approxi- 
mately 1.5% of their total weight during postcuring. For this reason, 
the following precautions were taken to control moisture absorption: 

• The specimens were dessicated immediately after postcuring. 

(The relative humidity in the dessicator was monitored and 
kept below 20%. ) 

• All specimens were tested within a month after postcuring when 
possible. 

• The specimens tested at room temperature were surrounded with 
dessicant that was kept in the oven space. 

In conclusion, our own moisture absorption experiments, together with 

a careful study of the literature available, let us believe that the 

precautions taken were sufficient to keep the moisture level well below 

.3% WT, a level that is believed to be negligible. 

5 . 3 Results and Discussion on 934 Resin 

Two 934-resin panels were fabricated by Graftech as detailed in 
Appendix A. These panels, labeled I and II, were rectangular plates 
of approximately 25" x 20" x .12", while one thick plate of 10" x 10" x 
.39" was made by Fiberite. The panels were cut into 8" x 1" 
rectangular specimens and labeled according to batch and position. In 
order to have as homogeneous a set of specimens as possible, each was 
subjected to a postcure cycle of 350°F (176°C) for 4 h as discussed in 
5.1. A screening experiment did not reveal substantial differences 
between the creep rates of specimens cut from different panels. 

The sequence of elevated temperature creep and creep-recovery 
tests, indicated in Fig. 5.12, was applied at a moderate stress level 
of a Q = 1600 psi (11 KPa). The inverse compliance at t = 1 min, i.e. 
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Time 

Fig. 5.12 Creep (1 hr) and Creep-Recovery (10 hr) Test Sequence at 
Different Temperatures. 

j Q /e (1 min) is plotted in Fig. 5.13. Moduli obtained from Instron 
tests are plotted on the same graph.* 

The creep and creep-recovery curves obtained during the test 
sequence given in Fig. 5.12 are very accurately modelled by the power 
law equation, 

D(T) = D Q (T) + D ] t n (2.19) 

and the result is shown in Fig. 5.14. 

*The Instron test results were obtained using a 2-in extensometer , 
whereas the creep and creep-recovery results were obtained using 
Micro-Measurements WK-06-1 25AD-350 strain gages. 
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Fig. 5.13 Neat Resin Modulus Versus Temperature. 


The temperature dependence of D Q (T), (Fig. 5.15), indicates a 
thermorheologically-complex material. Figure 5.15 also shows that the 
simple power law [Eq. (5.21)] breaks down above 300°F (149°C). The 
temperature dependence of D^(T) versus temperature is shown in Fig. 
5.16. The irregularity in this curve identifies the oj-transition. The 
transition is centered around 200°F (93.3°C), which is in agreement 
with the results obtained by Keenan et al . [79], i.e. 212°F (100°C), on 
this type of epoxy as well as those reported by Griffith [7] for the 
^-transition of the T300/934 composite system, 220°F (104°C). The 
value obtained by Yeow [5], 140°F (60°C), could be simply due to the 
combined effects of moisture absorption and uncertainties in the iden- 
tification of a small discontinuity in the thermal expansion versus 
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Fig. 5.16 Temperature Dependence of D^(T). 


temperature curve. 

The master curve, obtained by horizontal and vertical shifting 
of the data in Fig. 5.14 is shown in Fig. 5.17 while the amount of 
vertical versus horizontal shift is shown in Fig. 5.13. This latter 
figure can be approximated by two straight lines with slopes a and S 
respectively. The discontinuity in the slope is an indication of the 
change from one activation mechanism to another. It is well known that 
the horizontal shift in the glassy region can be represented by an 
Arhenius type of temperature dependence 


1 Ql 

°9 a j = 2.303 R 


I 1_ 

J ‘ Tr 


where aF is the activation energy (per mole), R is the universal gas 
constant, and T^ is an arbitrary reference temperature. Figure 5.19 
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Fig. 5.19 Horizontal Shift Versus the Reciprocal of Temperature. 

reveals a discontinuity in the log a^ versus T"^ curve. This again is 
an indication of a changing mechanism. The activation energy below 
220°C was found to be 4.76 kcal/gm-mole while above 220°F (104°C) a 
value of 23.50 kcal/gm-mole was found. 

Instron stress-strain curves (Fig. 5.20) revealed a brittle- 
ductile transition in the temperature range 200-300°F* (93-149°C). 

All curves have an initially linear stress-strain behavior. The 
values of stress and strain at the "linear-nonlinear transition point" 
are plotted versus temperature in Fig. 5.21. Based on these results, 
a test program was developed. At a given temperature, the same 
specimen was forced through a series of creep and creep recovery tests 

*Probably after the ^-transition has decayed. 
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Strain °/< 



Stress (Ksi) Stress (MPa) 





as indicated in Fig. 5.22. A fresh specimen was used for each new 
temperature. Test results obtained at five different temperatures, 
with at least eight subsequent creep and creep-recovery tests at each 
temperature, up till 320°F (160°C) did not show evidence of nonlinear 
viscoelastic behavior. As an example, Fig. 5.23 shows that the 
compliance versus time evolution at the highest temperature, 320°F 
(160°C), is stress independent. Similar results were obtained for the 
lower temperature tests. 

When we proceeded through the loading-unloading cycles (Fig. 5.22) 
we observed an accumulation of nonrecoverable strain (Fig. 5.24). The 
importance of this permanent strain accumulation process with respect 
to the viscoelastic deformation process can be judged when we replot 



Time 


Fig. 5.22 Test Sequence of Creep and Creep Recovery Tests at a Single 
Temperature and Different Stress Levels. 
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Log Time (min ) 

Fig. 5 .23 Comparison of Resin Compliance for Low (1.63 ksi (11.2 MPa)) 
and High (3.95 ksi (27.3 MPa)). 


Stress ( Ksi ) 


Fig. 5.24 Permanent Strain Accumulation Versus Stress (for 934 Resin 
at 246°F ( 11 9°C) ) . 










Fig. 5.24 with a vertical scale R, being defined as: 


total permanent set accumulated after k creep 

n _ and creep-recovery cycles lnnr 

the time dependent deformation at stress ° 

level 

This figure (5.25) shows that the permanent set accumulation tends to 
30% of the transient deformation at the highest stress levels. The 
nonlinearizing parameters g£ and a a became stress independent when we 
subtracted the accumulated strain (Figs. 5.26 and 5.27). This leads 
to the conclusion that the matrix material behaves linear viscoelastic 
at lower stress levels and linear viscoelastic-plastic at higher stress 
levels. 
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Fig. 5.25 Relative Importance of Permanent Deformation Compared to the 
Viscoelastic Deformation for Different Stress Levels. 
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Fig. 5.26 (a) Nonlinear Parameter g~ for 934 Resin with Subtraction 

of Accumulated Permanent Strain. 
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Fig. 5.26 (b) Nonlinear Parameter g 2 for 934 Resin Without Subtrac 

tion of Accumulated Permanent Strain. 
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Fig. 5.27 (b) Shift Factor for 934 Resin Without Subtraction of 

Accumulated Permanent Strain. 
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Chapter 6 

EXPERIMENTAL RESULTS AND DISCUSSION FOR T300/934 COMPOSITE 
6.1 Introduction 

Only creep measurements were made in our previous studies [4-13]. 
The Schapery nonlinear viscoelastic method described earlier could not 
be used as it requires both creep and creep recovery data. The present 
investigation was undertaken to obtain the proper information on 
T300/934 laminates such that the Schapery analysis could be developed. 
We also wished to compare the results of our measurements with new 
batches of T300/934 from a new source, Fiberite (see Appendix A for 
curing details), and earlier batches supplied by Lockheed-Sunnyvale . 

A further objective was to compare the nonlinear viscoelastic proper- 
ties of the T300/934 composite with those of the 934 resin given in 
Chapter 5. 

The plane stress constitutive properties of a linear orthotropic 
material is completely characterized by four independent properties. 
These may be the compliances, D^ , D^ = D 21 , ^>21 anc * ^66' ^he first 
two, D-|-| and D-^, can be determined from a tensile test of a uni- 
directional laminate tensile coupon with the fibers in the load direc- 
tion (see Fig. 6.1). D ^ and D^^ can be determined from tensile tests 
of a unidirectional laminate with the fibers transverse to the load 
axis and the fibers at 10° to the load axis respectively. Performing 
creep tests to express these parameters as a function of time permits 
the calculation of orthotropic linear viscoelastic response for an 
arbitrary direction of the fibers. 
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Fig. 6.1 Off-Axis Test Geometry and Internal Stress State. 

Normally, because the graphite fibers of a G/E laminate are very 
stiff and nearly time independent, creep tests to determine time 
dependent responses are only performed for the 90° and 10° cases 
respectively. 

If an orthotropic material is nonlinear, D ^2 is not only a func- 
tion of 02 but may be a function of a-j and as well [lb]. This gives 

rise to an interaction among stress components not present for a 
linearly viscoelastic orthotropic material. In order to account for 
this interaction effect, various methods have been employed usually 
based on the introduction of a flow rule. Such an approach is not 
always straightforward due to the possibility of various stress inter- 
actions at the microlevel. A simple method to account for interaction 
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effects is to replace the one-dimensional stress by the octahedral 
shear stress. This stems from the intuitive hypothesis that creep is 
caused by the shearing of the polymer molecules past one another. 


A somewhat more intricate approach was proposed for fiber rein- 
forced materials by Lou and Schapery [49] and was in turn used by 
Dillard [9], Their procedure utilized the octahedral shear stress only 
in the matrix phase as a nonlinearizing parameter. The relationship 


between the octahedral stress is the matrix t 


_m 

oct 


and the plane stress 


components (oi.c^jt-^) shown in Fig. 6.1 was given through a micro- 
mechanics model. The compliance model proposed by Dillard [9] was 


given as, 
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( 6 . 1 ) 


An alternate approach would be the use of the free energy expan- 
* 

sion of Eq. (3.18) which could accommodate higher orders of the strain 
or stress components or their invariants. Hahn and Tsai [83] used a 
fourth order strain expansion of the complementary energy density to 
handle shear nonlinearities. In other words, temperature was excluded 
from the Gibbs free energy of the system. Vie feel this approach should 
be extended and explored further. 

Nonlinear viscoelastic effect was not observed in the transverse 
direction for our T300/934 material. During the final compilation of 
the results it was realized that this result should be reexamined, 
since there is some indication that the nonlinear mechanism could 
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operate on such a long time scale that it becomes impossible to 
discriminate between materials response and experimental error on the 
basis of short time accelerated characterization results. 

6.2 Results and Discussion on [90]g s 

The laminate used in this investigation consisted of sixteen uni- 
directional plies, with about sixty percent volume of T300 graphite 
fibers in a matrix of Fiberite-934 epoxy resin. 

Tensile specimens were cut from a unidirectional plate by means 
of a diamond cutting wheel, postcured at 350°F (177°C) for 4 hours and 
dessicated as discussed earlier. A series of elevated temperature 
creep and creep-recovery tests were performed in the same manner as 
those for the resin shown in Fig. 5.12. The strain response was 
measured using back to back strain gauges to avoid the effects of bend- 
ing. Parenthetically, it is easily shown that an eccentricity (e) of 
1% of the depth (d) of the rectangular cross section as shown in Fig. 
6.2 causes a bending stress at the gauge location which is 6 % of the 
average stress . (This is discussed in more detail in Appendix C.) 



Fig. 6.2 Possible Bending of Tensile Specimen 
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The eccentricity in creep tension tests has been estimated by 
Isaksson [84] to be of the order of 5% of the depth (d) causing a 30% 
deviation of the average stress 

The change of the gauge factor with temperature was neglected 
because it was typical only 1.5% for a temperature change from room 
temperature up to 320°F (177°C). 

Another experimental necessity was to avoid heat leakage through 
the lower oven load train entrance as illustrated in Fig. 6.3. Efforts 
to insulate this entrance in the same way as the top load train 
entrance caused interference with the recovery data. The rod that was 
used to connect and disconnect the specimen with the lower and movable 
end of the load train needed to be completely free. A satisfactory 
solution for this problem was to place a heavy mass of metal at the 
lower entrance as shown in Fig. 6.3b. This "heat sink" contained an 



Fig. 6.3 Loading and Oven Insulation Details 
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opening trough which the connecting rod could move freely, but at the 
same time allowed a uniform oven temperature to be maintained. 

The transverse compliance as a function of time for different 
temperatures in the range 200-320°F (93-160°C) is shown in Fig. 6.4. 

An expanded vertical scale was used as compared with the similar resin- 
compliance plot of Fig. 5.3, because the continuous interruptions of 
the matrix by time independent filaments prohibited a large transverse 
creep strain buildup. 

Additional tests at different stress levels were performed for the 
transverse compliance as a function of temperature as shown in Fig. 6.4. 
However, these are not shown as the compliance appeared to be inde- 
pendent of stress level, i.e., was linear, for the time scale of our 
tests as will be discussed subsequently. 



Log Time (min) 

Fig. 6.4 Transverse or D 22 Compliance of T300/934 as a Function of 
Time and Temperature 
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The inverse compliance at t = 1 min is plotted versus temperature 
in Fig. 6.5. Similar results obtained by Yeow [5] and Griffith [7] are 
plotted on the same graph, together with a micromechanics prediction, 
based on the previously measured resin modulus as a function of 
temperature. 
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The micromechanics equation used to produce Fig. 6.5 was the 
Halpin-Tsai equation [54] 


E 2 (t,T) = E m (t,T) 


1 + S*n*v, 


1 - n*v. 


( 6 . 2 ) 


where 


E f /E m - 1 


in which 5 is a measure of reinforcement and depends on the fiber 
geometry, packing geometry and loading conditions as discussed in [54]. 
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Equation 6.2 is valid only for an elastic material but can be extended 
to the viscoelastic problem using the so-called "quasi-elastic" 
approach. It means that if the solutions to the elastic field equations 
are found, an approximate viscoelastic solution can be obtained by re- 
placing elastic constants by time dependent moduli. Sims [85] showed 
that for a creep time of several hours, this quasi -elastic approach 
deviates little from the more theoretically justified Laplace transform 
inversion approach. In Fig. 6.5 it was also tacitly assumed that 
E(t) = 1/D(t), which, according to Aklonis and Tobolsky [86], does not 
introduce an error larger than our experimental uncertainty in deter- 
mining the compliance. 

As may be observed, good agreement was obtained between the Hal pi n- 
Tsai micromechanics equation and current and previous experimental work. 
Such is the case even though the model is based on isotropic properties 
on both fibers and matrix even though the graphite fibers are known 
to be very anisotropic. This confirms a statement by Whitney [87] 
which indicates that a straightforward substitution of appropriate 
isotropic constants appears to be adequate to account for transversely 
isotropic fiber properties. 

Micromechanics solutions based upon an exact or approximate solu- 
tion to the field equations have been developed, however, which include 
complete anisotropy of the fibers. These include the upper and lower 
bound method of Bulavs et al . [90] and Hashin [91]. Experimental 
studies have been conducted by Ishikawa et al . [88] and Thompson [89] 
which show the Thornel 300 high strength fibers to be greatly aniso- 
tropic. These results are shown in Table 6.1. 
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Table 6.1 



Ishikawa et al . [88] 
(T-300A) 

Thompson [89] 
(T-300, wyp 15) 

E 11 

32.5 

(MSI) 

32 ± 3 (MSI) 

E 22 

2.93 

(MSI) 

/ 

v 12 

.2 


/ 

v 23 

.42 


/ 

G 1 2 

6.82 

(MSI) 

/ 


A comparison of Halpin-Tsai equations (6.2), Hashin's upper and 
lower bound method with anisotropic properties, Hashin's equations with 
equivalent longitudinal (E^ ) and transverse (Ef t ) fiber properties, 
Bulav's square array method and our current experiments are given in 
Fig. 6.6. As may be observed, Whitney's conclusion is once again con- 
firmed. 

Hashin's upper and lower bound results are shown in Fig. 6.7 as 
a function of the ratio of longitudinal to transverse fiber moduli. 

Equation 6.1 was modified to include the micromechanics approach 
given by Eq. 6.2. Using the quasi -elastic approach mentioned pre- 
viously predictions of the time dependent transverse modulus were made 
and these are given together with experimental results in Fig. 6.8. 

Good agreement was obtained which would indicate a simple time dependent 
micromechanics analysis to be a reasonable approach to represent our 
data over the time range involved. 
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Fig. 6.7 



Evolution of Upper and Lower Hashin Bound as a Function of 
the Ratio of Longitudinal and Transverse Moduli 
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Time ( min ) 


Fig. 6.8 Comparison Between Experimental (-) and Micromechanics (a) 
Prediction of E 22 Versus Time (T300/934 


An attempt was made to model the creep response with a Schapery 
type power law equation 


a (T).„(T) 

D(t) = g* T V + — — zr ^ — D. *t n 
v 0 n I 


(6.3) 


a 

where all quantities are as defined previously in Equation 3.49. The 
upper index T indicates that the nonlinearizing functions are evaluated 
in the transverse fiber direction. 

A series of 40 (one hour) creep and (ten hour) creep recovery tests 
at eight different stress levels and five different temperatures, using 
a fresh specimen at each new temperature, did not show any dependence of 
the transverse compliance on the stress level. In other words, 


(T) = „(T) , „(T) = 


= g 


= a = 1 


’0 J 1 g 2 a 

which means that the transverse response was linear within our time 
frame and within experimental error. 
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The stress-strain diagram (Fig. 6.9) obtained from an Instron 
test at 100°F (37°C) also shows little dissipation and a fracture strain 
of about 5000 microstrain. 



0.0 0.25 0.50 


Strain (%) 

Fig. 6.9 Stress-Strain Curve (Solid Line) for [90]g s T300/934 Laminate 
at Room Temperature (0.05 in/min) Compared with Linear 
Response (Dashed Line). 

Equation (6.3) reduces to 

D(t) = D q + D 1 t n (6.4) 

Figure 6.10 shows the temperature dependence of the first term (D Q ) in 
Eq. (6.4). It is clear that the dependence is almost perfectly linear 
in the temperature interval which we studied. The temperature dependence 
of the coefficient D-j has a pronounced nonlinear variation as can be 
seen in Fig. 6.11. Notice that the w-mechani sm which was identified in 
the neat resin (Fig. 5.2) is also clearly visible in Fig. 6.11. However 
the minimum in the curve is now approximately 20°F (* 10°C) higher. 
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Fig. 6.10 Instantaneous Transverse Compliance as a Function of 
Temperature 
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Fig. 6.11 D, (T) as a Function of Temperature 
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Such differences have been reported for the primary transition by a 
number of investigators [76,79,92], but the cause is unclear. 

The long time creep data plotted in Fig. 6.12 reveals that our 
experimental time scale for short time creep (1 hour) was probably too 
short to be able to discriminate between stress dependence of the 
compliance curve and small experimental irregularities such as mis- 
alignment in the load train, fluctuation in temperature, etc. Or, in 
other words, the shorter the test period, the more precise the test 
conditions need be in order to allow stress dependent discrimination of 
the compliance curve. Thus our earlier conclusion about the linear 
viscoelastic or stress independence of our transverse compliance data 
for T300/934 needs to be reexamined over a longer time scale. However 



Fig. 6.12 Comparison of Transverse Compliance Versus Time for Moderate 
(2.76 Ksi) 19 Mpa) and High (5.18 Ksi = 35.7 Mpa) Stresses 
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for a first approximation, linear would be a reasonable assumption. 

6.3 Results and Discussion on [1 0]ss 
6.3.1 Experimental Background 

An efficient intralaminar (in plane) shear characterization tech- 
nique known as the 10-deg off-axis test has been proposed and analyzed 
by C. C. Chamis and J. H. Sinclair [93]. Their analysis shows that the 
major contribution to fracture comes from intralaminar shear, i.e. the 
configuration activates a long continuous shear path, along which load 
diffusion into the reinforcing fibers occurs. One difficulty with the 
off-axis test is that shear coupling effects are a function of the load 
to fiber angle or the aspect ratio. The Pagano-Halpin [94] shear 
coupling analysis was used to generate Fig. 6.13. Their analysis sug- 
gest that for an aspect ratio or off-axis angle of 10°, the shear 
modulus, ^12 Approx ’ 1S a PP rox i ma tely 10% higher than the true shear 
modulus G^ 2 * This is confirmed by the inverse shear compliance versus 
temperature plot of Fig. 6.14 which indicates a measured value some- 
what above .8 MSI at room temperature, while a widely published value 
is closer to .7 MSI . 

Two back-to-back 45 deg rosette strain gages [WK-06-062R-350] 
were used to eliminate out of plane bending from test results. The 
orientation of the gages with respect to the fiber axis was checked 
under the microscope. Misalignments were taken into account during 
subsequent data reduction. Figure 6.15 shows the orientation of the 
strain gage rosette as well as the manner in which aluminum tabs were 
adhesively bonded to each specimen. 
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Fig. 6.13 Ratio of Actual to Apparent Shear Modulus in Off-Axis Test 
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Fig. 6.15 Specimen, Strain Gage and Tab Configuration. 

Each specimen was postcured with the same temperature cycle as 
were the transverse and neat resin coupons. 

The shear compliance D gg (t) for a 10° specimen loaded in tension 
as illustrated schematically in Fig. 6.1 is given by the following 
expression 

w.h. (-.598 e (t) + 1.879 e d ,(t) - 1.282 e (t)) 

D 66 (t) ’ " 

where e x (t), 645 (f) and 6 y^ are the strains measured a1on 9 the 
longitudinal , 45° and transverse directions of the specimen, respectively. 
The angle between the measured strain directions and the fiber-direction 
of the material are 10°, 45° and 100° respectively. W and h represent 
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the width and thickness of the tensile coupon. The strains measured 
along the x and 45° direction were positive, while the strain measured 
along the y-direction was negative. It should be noted, of course, that 
if the rectangular rosette is mounted with the 90° and 45° arms dif- 
ferent than shown in Fig. 6.15, different signs will be obtained and 
the preceding equation for D gg must be altered appropriately. The shear 
compliance D gg is thus easily obtained through simultaneous measurement 
of e x (t), e^ g (t) and e^(t) but caution should be taken with respect to 
the gage orientation. 

6.3.2 Experimental Results 

The shear compliance versus log t curves obtained at four dif- 
ferent temperatures are shown in Figs. 6.16a, b, c and d. They clearly 
indicate the behavior to be dependent upon stress. In the context of 
nonlinear viscoelasticity, as discussed earlier in Chapter 3, this means 
that a set of molecular processes which are sensitive to stress must 
also contribute to the observed deformation in addition to those nor- 
mally associated with linear viscoelasticity. 

A definite accumulation of nonrecoverable permanent strain was 
observed during the creep and creep-recovery test sequences. The amount 
obtained at 246°F (119°C) is shown as an example in Fig. 6.17. Thus the 
molecular structure moves out of the positions which it occupied on 
release of the stress but does not reach its former spatial arrangement 
under the action of the restoring forces. One possibility is that this 
former spatial arrangement is unattainable due to a "plasticity" type 
structural change in the matrix phase. Thus, for our material it seems 
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Fig. 6.17 Permanent Strain Accumulation During 10 Deg. Off Axis Test 
at 246°F (11 9°C) . 

not possible for nonlinear viscoelasticity to occur without involving 
other irreversible deformation mechanisms. This conclusion disagrees 
with results obtained by Ferry [28] on PMMA. Shear creep deformation 
up to 5 % were fully recoverable although the linear viscoelastic limit 
was below 1%. 

A number of investigators mechanically condition specimen prior 
to creep testing. The reason is that during the first creep test a 
certain amount of deformation takes place which does not recover upon 
removal of the stress. In subsequent creep tests on the same specimen 
at the same or smaller load there is no such nonrecoverable creep. The 
stress-strain relationship thus becomes repeatable. The material is 
therefore considered "mechanically conditioned." Leaderman [22] showed 
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that the permanent set in fibers disappears when the fiber was steamed 
under no load and then dried. Thus no real structural change or 
plasticity was introduced by the conditioning process. 

This author fears, however, that a significant volume of the 
matrix material might be yielded by the very inhomogeneous stress fields 
caused by the geometrical nature of the composite material. Evidence 
of such a matrix yielding phenomenon in shear deformation, based on a 
finite element analysis, was published by Foye [95]. His results shown 
in Fig. 6.18a give some indication of the progress of the plastic zone 
boundaries due to shear loading in repeating elements of a graphite- 
epoxy composite with 50% fiber by volume. Yielding initiates at the 
interface in an area enclosed by plastic zone boundary number 1. The 
plastic zone boundaries numbered 2 through 8 correspond to increasing 
postyielding load levels. Figure 6.18b, on the other hand, shows that 
in transverse deformation, yielding first occurs at a point in the 
matrix midway between two adjacent fibers. Foye also showed that initial 
yielding of the composite during the transverse deformation mode occurs 
at slightly lower average normal stress levels than the bulk matrix. 

In the shear deformation mode the composites yielded at much lower 
average shear stress levels than the bulk matrix. 

In conclusion, preconditioning was not made part of the testing 
program since it was feared that the stress-strain relationship would 
become repeatable at the expense of changing the matrix phase to an 
entirely new material in the sense of a plasticity type shakedown model. 
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Fig. 6.18a Inhomogeneous Stress Distribution Due to Shear Deformation 
Due to Shear Deformation After Foye [95]. 
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Fig. 6.18b Inhomogeneous Stress Distribution Due to Transverse Deforma- 
tion After Foye [95]. 

The creep response in the shear deformation mode was modeled with 
the power law 



200 


(6.5) 


D(t) = 


g(^D + 

y 0 0 


3 {S) Q (S) 

h g 2 


D ] -t' 


T 

The upper index (S) indicates that the nonlinearizing functions are 
evaluated in the shear mode. Plots of the nonlinearizing functions 
e gj^, 9^ > anc * a T obtained at 246°F as a function of stress are 
shown in Figs. 6.19a to 6.1 9d. Similar curves are obtained at higher 
temperatures. It can also be seen that gj S ^ and g^ vary in such a way 
with stress level, that their product is nearly unity and equation (6.5) 
for creep response becomes approximately 


oM - g' s) d„ + -i t" 


( 6 . 6 ) 


6.3.3 Time-Stress Superposition Principle (TSSP) 

A very useful viscoelastic concept is that time and stress are 
equivalent for describing viscoelastic behavior [7]. This is due to the 
fact that an increased tensile stress shifts the creep spectra towards 
shorter times as shown in Fig. 6.20a. Thus, as shown in Fig. 6.20b, 
creep observed for short times at a given stress x^ is identical with 
creep observed for longer times at a lower stress x Q provided that these 
data are shifted on a logarithmic time axis. Similarly, portions of the 

creep curve for x = x Q can be observed for stresses x = x^xg as 

indicated by the curved segments shown in Fig. 6.20b. These curved 
segments can then be shifted along the log time axis to construct a 
composite curve or the sigmoidal master curve applicable for the given 
stress level x , extending over many decades of time. 
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Fig. 6.19a Nonlinear Parameter g o (x). 
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Fig. 6.19b Nonlinear Parameter g-|(x). 
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Fig. 6.19c Nonlinear Parameter g 2 (x). 
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Fig. 6.1 9d Nonlinear Parameter a (x). 
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Log Time 


Fig. 6.20a Influence of Stress on Retardation Spectrum. 



Fig. 6.20b Dependence of Compliance D on log t and Stress t. 
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The logarithm of the shift factor a t for a particular curved 
segment is the horizontal shift necessary to allow it to join smoothly 
into the master curve. Each of the curved segments in Fig. 6.20b can 
be described by 


e 



D o + 


9l (S) 9< S > c 


(6.7) 


This same equation has a limited capability to represent a master 
curve using the form. 


D = 





C 



( 6 . 8 ) 


The reason for this limitation is due to the approximate nature of the 
power law used and was discussed earlier in Chapter 2. 

The time-stress-superposition principle (TSSP) described above 
was used to convert the data of Figs. 6.16a, b, c and d to master 
curves for each temperature and are given in Figs. 6.21a, b, c and d. 
Also, equation (6.8) was used to analytically determine each master 
curve and the resulting shift parameters. The symbols represent the 
master curves produced graphically while the solid line represents the 
analytically produced master curve. Comparison of the shift factors 
found by the two methods are shown in Figs. 6.22a, b and c. A three- 
dimensional version of these results is shown in Fig. 6.23. As may be 
observed, excellent agreement was obtained between graphically and 
analytically produced master curves. 
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Fig. 6.21a S gg Master Curve for T300/934 at 246°F (119°C). 
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Fig. 6.21b S gg Master Curve for T300/934 at 300°F (149°C). 
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Fig. 6.22a Comparison of Graphical and Analytical Obtained Horizontal 
Shift Factor .10 deg - off - axis for 246°F (11 9°C ) . 
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Fig. 6.22b Comparison of Graphical and Analytical Obtained Horizontal 
Shift Factor .10 deg - off - axis for 300°F (149°C). 







Fig. 6.22c Comparison of Graphical and Analytical Obtained Horizontal 
Shift Factor . 10 deg - off - axis for 335°F (168°C). 



Fig. 6.23 
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Graphical and Computed Shift Function Surface for T300/934. 
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Often, a vertical shift as well as a horizontal shift is neces- 
sary to produce a smooth master curve. In fact, even a rotation of 
individual compliance curves may be necessary. The type of shifting 
needed depends on the behavior of the limiting compliances D u and 
as a function of stress as shown schematically in Fig. 6.24. A 
vertical shift of D x - D U T ° produces the dotted durve D x (t) which can 
be superimposed on D x °(t). 

A rotation would be necessary when the limiting compliances are 
affected differently by the stress intensity, i.e., when 

C ~ D T 

— — f 1 (6.9) 

D x ° - D T ° 

R u 



Stress Log Time (min) 

Fig. 6.24 Illustration of the Stress-Shifting-Procedure (Based on 
McCrum et al . [96]) . 

In our case, a small vertical shift was found necessary at 246°F. This 
shift however was nonexistent at higher temperatures as shown in Figs. 
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6.25a, b and c, where the vertical shift is plotted versus the hori- 
zontal shift. 

It was thought initially that the true vertical shift 
" C°), given in Fig. 6.24, was equivalent to (g Q - 1)D Q . This 
however leads to an overprediction due to the fact that the sigmoidal 
compliance curves change location with increasing stress. The impli- 
cation is that even when the limiting compliance values are completely 
insensitive to stress (which turns out to be the case in this investi- 
gation) we still would record a g Q value different from unity since the 
retardation spectrum involves rapid creep on the time scale of micro- 
seconds. In conclusion, very accurate short time limiting compliance 
data is needed to obtain a proper estimate of the true vertical shift. 

As a final conclusion, a long time control test remains the only 
criterion for checking a particular accelerated test method. These 
long term creep tests are currently being carried out and will dictate 
the adaptations necessary to achieve a proper mathematical modeling for 
the prediction of long term properties from short term test results. 
These final modeling equations should in some way reflect the mechanism 
of inelastic phenomena operating at the microstructural level. 
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6.25a Cross-plot of Horizontal and Vertica 
Shift at 246 Deg F (119°C) 
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6.25b Cross-plot of Horizontal and Vertica 
Shift at 300 Deg F (149°C) 
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Fig. 6.25c Cross-plot of Horizontal and Vertical 
Shift at 320 Deg F (160°C) 
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Chapter 7 

EXPERIMENTAL RESULTS AND DISCUSSION ON DELAYED FAILURE 
AND LAMINATE RESPONSE 

7.1 Introduction 

The usual concept of elastic deviatoric strain energy in which 
the Tsai -Hi 11 anisotropic failure criterion is based cannot be carried 
over directly to viscoelastic materials without modifications to 
account for dissipation mechanisms. However, a better method to in- 
clude dissipation as well as the entire history of deformation is the 
concept of free energy which can be used to define a "critical" stress 
state at the point of failure. A delayed failure model based upon 
these concepts was developed in Chapter 4. 

The delayed failure model requires information on both creep 
and on creep rupture. An existing creep rupture data base for [903g s 
and [60]g s orientations of T300/934 was available from the work of 
Griffith [7]. [90]g s creep data is given in Chapter 6. The needed 

[60]g s creep response to compliment the [60] gs rupture data can be 
inferred from the [90]g s and [10]g s data given in Chapter 6. However, 
it was decided to obtain the needed [60]g $ creep response from new 
specimen because inferred response would introduce a hypothetical 
octahedral stress interaction factor, t™ t . 

7.2 Delayed Failure 

7.2.1 Analysis on [90]g $ Results 

The compliance versus time curve at 320°F (160°C) and 1580 psi 
(10.88 MPa) (Fig. 6.8) was fitted with a six-element Kelvin model. 
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using a collocation scheme proposed by Schapery [97]. The resulting 
compliance D(t) can be expressed as: 

D(t) = .71 + .041 [1 - exp (-t/.Ol)] + 

+ .041 [1 - exp ( -t/1 ) ] + .060 [1 - exp ( -t/1 0 ) ] 

+ .049 [1 - exp (-t/1 00)] + .164 [1 - exp (-t/1000)] 

+ .158 [1 - exp (-t/10000)] x 10" 6 psi" 1 (7.1) 

Based on this identification scheme, we can use Eqs. (4.9 - 4.11) 
to calculate the total creep energy, stored free energy and dissipated 
energy as a function of time as shown in Fig. 7.1. A gradual increase 
of the total dissipation with time originates from the energy con- 
sumption in the subsequent viscous elements. This dissipation has to 
be subtracted from the total creep energy, which is in this case equal 
to the potential energy loss of the external loading, in order to ob- 
tain the total stored free energy. 

An identical evolution of the free energy was also obtained based 
on the power law parameters (D Q ,D^,n) which made it possible to 
evaluate Eq. (4.17) as a function of time. This procedure has a 
definite advantage in that it avoids a collocation scheme. 

A comparison of the experimental creep-rupture results and the 
free energy based prediction shows that very good agreement was ob- 
tained as shown in Fig. 7.2. The full line (-•“) represents a given 
constant free energy line. Thus a creep rupture result is obtained 
from creep compliance data at the expense of only one additional 
parameter (w). Earlier attempts, in which we tried to correlate the 
total creep-energy with creep to rupture data, gave very poor agreement. 
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Fig. 7.1 T300/934: Evolution of Total Creep Energy, Total Stored 

Energy and Dissipated Energy as a Function of Time for 
[90 ] 8s at 320 deg F (160° C). 



Log Time (min) 


Fig. 7.2 T300/934: Comparison of Free Energy Based Creep to Rupture 

Prediction with Experimental Results Obtained on [90] Qc 
at 320 deg F. 8S 
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7.2.2 Analysis of [603g s Results 

Long time creep and creep recovery data were obtained on a 60-deg 
off axis specimen. The details on postcuring, strain gaging and test- 
ing were identical with those already discussed in Chapters 5 and 6. 

A plot of the compliance change with log t at a stress level of 
72% of ultimate is shown in Fig. 7.3. A result obtained by Griffith 
[7] on the same orientation but at 55% of ultimate is also shown in 
this figure. The curvature of these plots is clearly affected by the 
stress intensity. In addition, however, a rather large vertical shift 
would be needed to superimpose both curves. This could be due to 
matrix plasticity, since a rather large permanent deformation remained 
after creep recovery (Fig. 7.4). Martirosyan [98] gave evidence of 
such results which are reproduced in Fig. 7.5. The nonrecoverable 
creep deformation is seen to behave as an increasing-decreasing func- 
tion with fiber angle. These bell -shaped curves are parameterized by 
a number that indicates the stress intensity as a fraction of ultimate. 
Obviously a large nonrecovered strain was observed which varies very 
nonlinearly with stress intensity. 

The creep response curve at 72% of ultimate was again modeled 
with a six element Kelvin model, with an additional free dashpot in 
series, as shown in Fig. 7.6. 

The free energy at rupture is seen to be constant for delayed 
failures in excess of 30 min. as shown in Fig. 7.7. Exactly the same 
free energy evolution with time was found again through the power law 
approximation. The power law parameters were evaluated, however, 
after subtraction of the free dashpot contribution. 
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Long Time Creep Compliance of [60] o „ T300/934 at 320°F 
(160°C). 8S 
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Fig. 7.4 Typical Creep Curve for [60]o c T300/934 at 320°F (160°C). 
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Fig. 7.5 Nonrecoverable Creep of Glasfiber Reinforced Plastic with 
Respect to Fiber Angle and Percent of Ultimate Stress. 



Fig.. 7. 6 Kelvin Element Used to Model Creep Response 
for [60]g s at 320°F (160°C). 
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Fig. 7.7 Free Energy at Rupture for [60] Qc T300/934 
at 320° F (160°C) . as 

A reasonable free energy based creep-rupture prediction was 
obtained, as seen in Fig. 7.8. Final conclusions however cannot be 
made yet, due to the limited number of experimental data points which 
are currently available. 

7.3 Creep Response on [90/±45/90]2 S * Compared with Creep response of 
[90] 8s and [±45] 4s at 320°F (160°C) 

It can be expected, that a [90/±45/9032 S laminate configuration 
will operate at a lower creep rate than a [±45] 4s , since the latter 
contains only two fiber orientations which allows scissoring action 
as may be visualized in Fig. 7.9. The former contains three fiber 
orientations and thus a fibre truss network is formed as shown in 
Fig. 7.9 which prevents scissoring. 
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7.9 Schematic of [±45] and [90/±45/90] Laminates with Fiber 
Truss Network in Latter. 
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Figure 7.10 compares the creep responses of three laminates. It 
is qualitatively clear that the [±45] 4$ fiber orientation shows the 
most transient creep. This can also be quantified when we compare the 
power law coefficients, D-j , which are a function of stress, as given 
in Figs. 7.11. A comparison of the [90/±45/90 ] 2$ laminate results 
with those for the [90] gs and [±45] 4s laminates reveals that the creep 
rate for the former is an order of magnitude less than that for the 
other orientations. The initial compliance however decreased only by 
35% as compared to the [±45] 4$ orientation. 



0.0 I .0 2.0 

Log Time ( min) 

Fig. 7.10 Creep Compliance for Three T300/934 Laminates. 
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Three Different Laminates at 320°F (160°C) 


122 


Nonrecoverable Deformation 
( microstrain ) 





In conclusion, adding 90' s to a ± 45 orientation increases the 
modulus by only 35% but decreases creep sensitivity by an order of 
magnitude. 

Prior to testing the [±45]^ s laminate, we also put a strain gage 
transverse to the loading direction. This information allows us to 
obtain the shear compliance [99] which is given by 
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2(6 x - V 


(7.2) 


It was observed, however, that e (t) and e (t) had approximately the 

x y 

same magnitude. This led to untolerable large errors and even in- 
consistencies in the creep recovery data. 


7.4 Numerical Prediction of [90/±45/90] 2s Laminate Behavior 

The incremental nonlinear viscoelastic lamination theory-program, 
which was written and documented by Dillard [7 ] during an earlier 
phase of our research program, was used to generate the creep response 
of a [90/±45/90]2 $ laminate. 

A comparison of the numerical prediction with the experimental 
result at 70% of ultimate at 320°F (160°C) is shown, on an expanded 
scale, in Fig. 7.12. Reasonable agreement between predicted and 
measured compliance was obtained. Creep rupture-time, however, was 
even more seriously underpredicted, indicating a basic deficiency of 
the numerical approach which lies in its inability to reflect the 
microstructural effects of the laminate. This microstructural behavior 
was deduced using edge replicas [100], as shown in Fig. 7.13. An 
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Fig. 7.12 Creep Compliance of a [90/±45/90]« T300/934 Laminate at 
320°F (160°C) . “ 
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Edge replica after 15,840 minutes 
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Fig. 7.13 Edge Replicas of [90/±45/90] 2s Laminate. 
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attempt to model the stress transfer mechanism involved by means of 
a shear-lag model [101], which was simple enough for subsequent 
implementation in the laminate analysis program, turned out to be un- 
successful . 

Comparison of the creep response at 97% of ultimate, also as 
shown in Fig. 7.12, with the response at 70% of ultimate indicates a 
large vertical shift. The edge replica, obtained after unloading, 
however, revealed a very different crack spacing between results at 
the two stress levels. The crack spacing in the 90-deg. oriented 
lamina is considerably less than in the previous case and some 
cracking in the 45-deg. oriented lamina became visible at the higher 
load level. 

Figure 7.14 shows the irregularities in the strain reading at 
the beginning of the creep rupture experiment at 97% of ultimate. 

This indicates a sequence of events at the microstructural level, 
which could be deterministic in the selection process of the stress 
transfer mechanism. 

More understanding in this area should be generated, including 
the answer on the question of whether or not crack healing sets in 
under certain unloading or strain-recovery conditions. 
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Chapter 8 

SUMMARY AND RECOMMENDATIONS FOR FURTHER STUDY 

Different constitutive formulations able to model material non- 
linearities in viscoelastic systems have been discussed. Nonlinear 
viscoelasticity models are often quite complicated, creating great 
experimental difficulties for the determination of material parameters. 
It was concluded that the approach advocated by Schapery is sufficiently 
general to account for stress-activated nonlinear behavior in polymers 
and composites. Particularly the following experimentally observed 
phenomena were easily modeled: 

• An initial elastic creep-recovery which is greater than the 
initial elastic creep response 

• A recovery behavior which is overpredicted by linear visco- 
elastic considerations (i.e. the actual recovery curve in the 
nonlinear range is flatter than in the linear range) 

Creep and creep-recovery experiments which were carried out on 
the neat 934-resin did not reveal any stress dependent change in 
materials response other than plasticity. The analysis of time- 
dependent response, which is exclusively connected with specifics of 
the molecular structure revealed the existence of a secondary transi- 
tion-temperature located between 200 and 250°F ( 93-1 21 °C ) . 

Essentially the same series of experiments were carried out on 
the T300/934 graphite-epoxy system. The transverse compliance D 22 (t) 
was found to be very stress insensitive. The shear compliance Dgg(t), 
on the other hand, turned out to be a sensitive function of stress. 
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The role of yielding and subyielding of the matrix phase, due to dif- 
ferences in loading path (transverse loading vs. shear loading) was 
discussed. Two uncoupled equations were fit to the data D^U) anc * 

Dgg(t) respectively. The calibration of the material parameters 
(D o ,C,n,g o ,g^ .g^ja ) was done by a computer routine. It is well known 
however that mathematical models fit data well when used on the stress 
paths which were applied for parameter calibration purposes. In order 
to account for stress interaction along general loading paths, the 
principle of octahedral stress, was borrowed from nonlinear approaches 
in isotropic viscoelasticity. The validity of this approach however 
needs to be carefully checked. The power law creep exponent (n), was 
found to be .23 and independent of stress and temperature. The 
secondary transition, which was identified in the resin, was also 
identified in the composite, but at a slightly higher temperature. 

Short time data could be shifted graphically into a unique "master 
curve." The horizontal shift was found to be in very good agreement 
with the analytical shift (log a ). The vertical shift was very small, 
and became even nonexistent at higher temperatures. It was shown that 
vertical shifting is related to the behavior of the limiting compliances. 
An interesting point in this context is that Urzhumtsev [102] presents 
evidence that the time-stress-superposition principle, when used to 
shift short term compliance curves obtained on porous polyurethane is 
in excellent agreement with a long time control experiment whereas the 
time-temperature superposition leads to very serious overpredictions 
of long time compliance data. 
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It should be noted that no mechanical conditioning was applied, 
since conditioning should be used to isolate delayed elastic response, 
which is something that has never been proven for fiber reinforced 
composites (at least not to this author's knowledge). 

Creep compliance data were merged into a free energy-based 
failure criterion and reasonable creep-rupture predictions were ob- 
tained. The creep response of a [90/±45/90]2 S laminate was found to 
be an order of magnitude less than the creep response of [90]g s and 
[±45] 4$ , while the modulus increased only by a factor of 1.35. A 
reasonable creep compliance versus time prediction was obtained by means 
of a nonlinear viscoelastic laminate analysis program. It is felt, 
however, that a damage model that operates on the minimechanics level 
of the problem would give a better agreement than the cumulative 
damage model which was actually used [9]. Promising ideas in this 
direction were published recently by Nuismer et al . [103]. They 
developed a realistic and simple scheme that accounts for influences 
on the constitutive behavior of given matrix dominated laminas through 
the lamina-laminate interaction after damage initiation. Another, more 
fundamental, approach that describes the influence of damage through a 
continuum vector field has been pursued, within a continuum thermo- 
dynamics framework by Talreja [104]. 

Based on this current research effort, the following recommenda- 
tions for further study seem appropriate: 

• The feasibility of mechanical conditioning, in case of 
composites, should be thoroughly checked. 
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Test programs for boron-epoxy by Cole and Pipes [14] and for 
graphite polymide by Pindera and Herakovich [15] produced evi- 
dence that the octahedral shear stress is not a valid 
representation of stress mode interaction. Pindera and 
Herakovich showed that, in order to get access to the actual 
stress interaction, 3 off-axis tests were needed. An 
identification of the actual representation for stress mode 
interaction in case of graphite-epoxy T300/934 would at least 
give indications on its relative importance with respect to 
nonlinearities in transverse deformation and shear. 

The free energy based failure criterion together with the 
Tsai-Hill criterion, could be developed to predict delayed 
off-axis failure, based on transverse and shear compliances 
only. 

The necessity of transient temperature tests, in order to 
isolate the exact vertical shift, and thus obtain maximum 
predictive power, should be explored. 

Reversibility is always dependent on the choice of a well- 
defined easily reproducible reference state. In composites, 
however, initial stresses arise as a result of physio- 
chemical processes, technological operations and mechanical 
and thermal incompatibility. The natural state, i.e. a body 
free of stresses and strains, is nonexistent. The importance 
of an "eigenstress"-field on the viscoelastic behavior should 
be checked. 
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APPENDIX A 


Cure Procedure for T300/934 Composite Panels 


1. Vacuum bag - insert lay-up into autoclave at room temperature. 

2. Apply full vacuum. (A minimum of 25 inches of Hg.) Hold for 30 
minutes at room temperature. 

3. Maintain full vacuum throughout entire cycle. 

4. Raise temperature to 250 F (+5. - 10 F) at 2 ± 5 F per minute. 

*5. Hold at 250 F (+5. - 10 F) for 15 ± 5 minutes, apply 100 psi 

(+5. - 0 psi ) . 

*6. Hold at 250 F (+5. - 10 F) and 100 psi (+5. - 0 psi) for 45 ± 5 
minutes . 

*7. Raise temperature to 350 F (+10. - 0 F) at 2 ± 5 F per minute. 

*8. Hold at 350 F (+10. -OF) for two hours ± 15 minutes. 

*9. Cool under pressure and vacuum to below 175 F. 


*Pressure was 85 psi . 
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Curing Procedure for Resin Panel #1 


1. Resin heated until fluid and poured into heated mold. 

2. Mold vacuum bagged and inserted into autoclave. 

3. Full vacuum applied and held 30 min. 

4. Raise temperature to 250 F (+5. -10 F) at 2 ± 5 F per minute. 

5. Hold at 250 F (+5 -10 F) for 15 ± 5 minutes. Apply 100 psi 

(+5 - 0) 

6. Hold at 250 F (+5 -10 F) and 100 psi (+5 - 0 psi) for 45 ± 5 min 

7. Raise temperatuer to 350 F (+10 -0 F) at 2 ± 5 F per minute. 

8. Hold at 350 F (+10 -0 F) for two hours ± 15 minutes. 

9. Cool under pressure (85 psi) and vacuum to below 175 F. 


Fabrication Procedure for Resin Panel #2 

1. Mold heated to 200 F. 

2. Resin heated to 160 F and degassed in bell jar for 5 minutes. 

3. Degassed resin poured into one corner of tilted and heated mold. 
Air gun used to keep mold warm. 

4. Resin then cured in oven with same schedule as prepreg tape, i .e 
steps 4-9 above. 
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The Curing Procedure for Resin Panel #3 


934 Resin Cure Procedure - Solution Form 


1. Break off chunks of solid (cold) resin. 

2. Place in beaker of acetone or methyl ethyl keytone. 

3. Slowly warm and stir solution up to 80°F (26°C) and until all 
resin is dissolved to taffy consistency. 

4. Pour into teflon coated bread pan and place in vacuum oven. 

Heat to 150°F (66°C) under full vacuum. Observe until bubbles 
cease (about 2 hours). 

5. Remove from vacuum oven and put in 21 2° F (100°C) oven over night 
(roughly 3:00 p.m. - 7:00 a.m.). 

6. Remove from oven in morning. 

7. Put in oven at 325°F - 350°F (163°C - 177°C) for two (2) hours. 

8. Turn off oven and allow to cool gradually. 
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APPENDIX B 


with 

and 


G = F - Q m q m 


F - F <VV T > 


Q = — 

m 3 < V, 


Thus 


G = G(Q m ,q r ,T) 
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APPENDIX C 


INFLUENCE OF EXCENTRIC LOAD APPLICATION [105] 

Suppose that the load F is applied with a small excentricity e. 
This causes not only the homogeneous stress state which we intended 
to produce, but also the bending moment. F-e as is shown in Fig. C.l. 


F F 



Fig. C.l Superposition of Normal Stress and Bending Due to Excentric 
Load Application 


The stress, at the gage location a QL can then be written as, 


. - F + F-e 

a GL ' w-t u t 2 
w • t 

6 


a u [’ + 6 't 


with a u = > the uniform stress, which we intended to apply. 
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